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Abstract 

Refining previous work in |41l 1341 1371 1161 I15| . we derive sliarp pointwise bounds on behavior of 
perturbed viscous shock profiles for large-amplitude Lax or overcompressive type shocks and physical 
viscosity. These extend well-known results of Liu |27| obtained by somewhat different techniques for 
small-amplitude Lax type shocks and artificial viscosity, completing a program set out in |42| . As 
pointed out in |2til 127) . the key to obtaining sharp bounds is to take account of cancellation associated 
with the property that the solution decays faster along characteristic than in other directions. Thus, we 
must here estimate characteristic derivatives for the entire nonlinear perturbation, rather than judicially 
chosen parts as in 1371 115| . a requirement that greatly complicates the analysis. 

1 Introduction 

In the landmark paper [27], Liu established sharp pointwise bounds on the asymptotic behavior of a perturbed 
viscous shock profile, for small-amplitude. Lax type shocks and artificial viscosity. A long-standing program 
of the authors, set out in |42j . has been to extend to large-amplitudes, general type profiles, and physical 
(partially parabolic) viscosities results obtained by Liu and others for small amplitude Lax type profiles 
and artificial viscosity. For various results in this direction, see, e.g., In this paper, we 

achieve a definitive result recovering the full bounds of |27| for general Lax or overcompressive type shocks 
and physical or artificial viscosity, thus completing the program of jl^l- The analysis involves an interesting 
blend of the spectral techniques introduced by the authors to deal with large amplitudes, described in 
Sections 12 and 121 a delicate type of cancellation estimate introduced by Liu in the small- amplitude context, 
described in Section^ and sharp convolution estimates of the type developed in ^l^l^^^lE], taking 
account of transversality but not cancellation of two interacting signals. 

Consider a (possibly) large-amplitude viscous shock profile, or traveling-wave solution 

u{x — st); lim u{x) — u±, (1-1) 

X — >±QO 

of systems of partially or fully parabolic conservation laws 

ut + F{u), = {B{u)u^),, (1.2) 

a; e R, M, e M", B e M"^", where 

-(;'). (-) 

E M"^^, u^^ EW,r some positive integer, possibly n (full regularization) , and 

Re cr(b2) > 0>O. 
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Here and elsewhere, a denotes spectrum of a matrix or other linear operator. Working in a coordinate system 
moving along with the shock, we may without loss of generality consider a standing profile s = 0. 

Following 40 , we assume that, by some invertible change of coordinates u — > followed if necessary 

by multiplication on the left by a nonsingular matrix function S'(w), equations H1.2() may be written in the 
quasilinear, partially symmetric hyperbolic-parabolic form 

A°wt + Aw^ = {Bw^)x +G, (^^jj^ , (1.4) 

e M"^'*", w^^ e W, .T e M, t e M+, where, defining w± := w{u±): 
(Al) A{w±), All, A" are symmetric, A° > 0. 

(A2) Dissipativity: no eigenvector of dF{u±) lies in the kernel of B{u±). (Equivalently, no eigenvector of 
i(i")-i(u;±) lies in the kernel of B{A^)-^{w±).) 

(A3) S = ( H, G = (^), with Reb{w) > for some 6* > 0, for all W, and g{w.^,w^) = 0{\w^\'^). 



V' \9/ 

Here, the coefhcients of (|1.4|l may be expressed in terms of the original equation 1)1. 2|) . the coordinate change 
u — > w{u), and the approximate symmetrizer S{w), as 

i° S{w){duldw), A := S{w)dF{u{w)){du/ dw), 

B -.^ S{w)B{u{w)){du/dw), G = -{dSw^)B{u{w)){du/dw)w^. ^^'^^ 

Alternatively, we assume, simply, 

(Bl) Strict parabolicity: n — r, or, equivalently, 5R(t(_B) > 0. 

Along with the above structural assumptions, we make the technical hypotheses: 
(HO) F, B,w, S e Ci°. 

(HI) The eigenvalues of A* := Aii{Aii)~'^ are (i) distinct from 0; (ii) of common sign; and (iii) of constant 
multiplicity with respect to u. 

(H2) The eigenvalues of dF{u±) are real, distinct, and nonzero. 

(H3) Nearby u, the set of all solutions of Hl.l(l - H1.2|l connecting the same values u± forms a smooth manifold 

Remark 1. Structural assumptions (Al)-(A3) [alt. (Bl)] and technical hypotheses (H0)-(H2) admit such 
physical systems as the compressible Navier-Stokes equations, the equations of magnetohydrodymics, and 
Slemrod's model for van der Waal gas dynamics [Ml- Moreover, existence of waves u satisfying (H3) has 
been established in each of these cases. 




Definition 1. An ideal shock 

u{x,t) = r- (1.6) 

I W+ X > st, 

is classified as undercompressive, Lax, or overcompressive type according as i — n is less than, equal to, or 
greater than 1, where i, denoting the sum of the dimensions i_ and i+ of the center-unstable subspace of 
df{u^) and the center-stable subspace of df{u+), represents the total number of characteristics incoming to 
the shock. 

A viscous profile (11.11) is classified as pure undercompressive type if the associated ideal shock is un- 
dercompressive and i — 1, pure Lax type if the corresponding ideal shock is Lax type and i — i — n, and 
pure overcompressive type if the corresponding ideal shock is overcompressive and £ — i — n, £ as in (H3). 
Otherwise it is classified as mixed under-overcompressive type; see 
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Pure Lax type profiles are the most common type, and the only type arising in standard gas dynamics, 
while pure over- and undercompressive type profiles arise in magnetohydrodynamics (MHD) and phase- 
transitional models. In this paper, we restrict to the case of pure Lax or overcompressive type shocks. 

Finally, we assume that the profile satisfies a linearized stability criterion based on the Evans function. 
As described, e.g., in |T1 ^21 ESI ED ISl El E3 . the Evans function D{X), a Wronskian constructed from 
solutions of the associated eigenvalue equation, serves as a characteristic function for the linear operator L 
that arises upon linearization of H1.2|l about u{x). More precisely, away from essential spectrum, zeros of 
the Evans function correspond in location and multiplicity with eigenvalues of L 02] ■ It was shown 
in 1321 and j^, respectively for the strictly parabolic and real viscosity cases, that DLp ^ LP linearized 
orbital stability of the profile, p > 1, is equivalent to the Evans function condition 

(V) There exist precisely £ zeroes of D{-) in the nonstable half-plane Re A > 0, necessarily at A = 0, 
where £ as in (H3) is the dimension of the manifold connecting u_ and u+. 

Under assumptions (A0)~(A3) [alt. (Bl)] and (H0)-(H3), condition (V) is equivalent to (i) strong spectral 
stability, (j{L) C {Re A < 0}U{0}, (ii) hyperbolic stability the associated ideal shock, and (iii) transversality 
of M as a solution of the connection problem in the associated traveling- wave ODE, where hyperbolic stability 
is defined for Lax and undercompressive shocks by the Lopatinski condition of |2H1 EHl 1^111 [Z| and for Lax 
and overcompressive shocks by the analogous long- wave stability condition (Pii), below; see 021 E3 ESI ED] 
for further explanation. 

Remark 2. Stability criterion {V) has been shown to hold for general small amplitude Lax shocks |51 llll[THl 
1231 124L 1551 136| . and for large-amplitude shocks in such cases as Lax type waves arising in iscntropic Navier- 
Stokes equations for the gamma-law gas as 7 — > 1 ^5 , and undercompressive shocks arising in Slemrod's 
model for van der Waal gas dynamics (SHI (see [SH] for Slemrod's model). On the other hand, it has been 
shown to fail for certain large-amplitude and or nonclassical type shocks |1(J[ O ESI ■ More generally, 
condition (2?) can be readily checked by numerical calculation ^ 4^ El El 1171 E| ■ 



Setting A± := df{u±), T± cPf{u±), and B± := B{u±), denote by 

< a2 < ■ ■ ■ < a^i and < < • • • < (1-7) 

the eigenvalues of and A^, and l^ , r^ left and right eigenvectors associated with each a^, normalized 

so that {ljrk)± — Sj., where Sj. is the Kronecker delta function, returning 1 for j = k and for j ^ k. Define 
scalar diffusion coefficients 

[if (;JSr,)± (1.8) 

and scalar coupling coefficients 

7± :=(Zjr(r„r,))±. (1.9) 

Under this notation, hyperbolic stability (a consequence of the assumed (2?)) of a Lax or overcompressive 
shock profile u is the condition: 

(2?ii) The set {rf : af ^ 0} U {J^^ ^dx] i = 1, • . . , £} forms a basis for M", with ^dx com- 
puted at (5 = 0. 



Remark 3. For Lax profiles, {T>\\) reduces to the Liu-Majda condition 

det ({r± : af ^ 0} (u+ - u_)) ^ 0. 

Following 1251 127j . define for a given mass raj the scalar diffusion waves (pj {x,t;m~) as (self-similar) 
solutions of the Burgers equations 

vj,t + «7^j> - i^IvJ,xx = -^i((vjf)x (1-10) 

with point-source initial data 

Vj{x,-l) = mj6Q{x), (1.11) 



3 



Howard, Raoofi, and Zumbrun, Liu Estimates 



and similarly for ip'j'ix, t; m'j'). Given a collection of masses prescribed on outgoing characteristic modes 



< and aj > 0, define 



(1-12) 



a, <0 



aj>0 



In the setting described above, we will determine estimates on perturbed viscous shock profiles in terms 
of (j) and a refined collection of template functions (terminology following ,42j|; notation following [Liu97]). 

Definition 2 (Template functions.). Let 



(1.13) 



M^,t)^{i + \x\+t)-'{i + \x\)-'x 

^4(x,t)-(l + N+t)-V4^, 

where x denotes an indicator function on x ^ [a^t^a^t], and also 

^{^{x,t)^{l + \x-aft\+t^/^)-^'^ 
(x, t) = -01 (a;, i) - 

The goal of our analysis is to establish the following sharp pointwise description of asymptotic behavior, 
generalizing bounds obtained by Liu PY) for small-amplitude Lax type profiles with artificial viscosity B = I. 

Theorem 1. Assume (Al)-(A3) [alt. (Bl)], (H0)-(H3) and [V) hold, and u is a pure Lax or overcom- 
pressive shock profile. Assume also that u solves with initial data Uq and that, for initial perturbation 

uq uq — u, we have \uq\h5 < Eq and \uo{x)\, \dxUo{x)\, and |9^uo(x)| < Eq{1 + \x\)^^ , for Eq sufficiently 
small. Then, the solution u continues globally in time, with u — u G L^ (1 . Moreover, there exist a 
choice of \m^\, |(5*| < Ci?o and a function 6(t) (determined, respectively, by 1^3. 1\ and ^3.1S}) ). such that, for 

V 



Hx,t)\ < CEo[tJji+4'2 

\dMx, t)\ < CEo [t~^/H^i + ^2) +i'3 + ^i 
m)\<CEo{l+tri 



ix,t) 
^2) 



{x,t) 



(1.14) 



\m\<cEo{\+tr^ 

and, for all ^ 0, 

\{dt + afd:,)v{x,t)\ < CEo[t'\l+t)^/^iji'^ +t- 

for some constant C (independent of x,t and Eq). 
Remark 4. By Taylor's Theorem, 



■1/2 ^A^i 



(V'l +^2) +V'3 +V'4 {X,t), 



(1.15) 



.k(-S,+5(t) 



-,s. 



Is 



5{t)) = 0{m\''e-^\^\), 



(1.16) 



for fc = [resp. fc = 1], since ^j- and its derivatives decay at exponential rate (see Lemma Since the 
righthand side of p.l6|l for fc = [resp. fc = 1] is smaller than the right hand side of (|1.14|l fi) [resp. (ii)], we 
may conclude that 

decays at the same rate as v, i.e., u is well approximated by -;2''*+'^(*) + 0: a dynamically changing nearby 
profile with the same endstates u±, superposed with outgoing diffusion waves. This is a slight refinement of 



the asymptotic description u • 



introduced by Liu [23 EH- 
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Remark 5. In contrast to the modulus bounds given here, the estimates of j^T] are stated in each separate 
characteristic field, giving the further information that along characteristic rays the perturbation v lies mainly 
in the associated characteristic direction. The difference is rather subtle, however, and so we have opted for 
simplicity of exposition to omit this level of detail. The additional information can easily be recovered at 
the expense of additional bookkeeping; see Remark [T^ for further discussion. 

Liu's analysis involved essentially two main ingredients, which were strongly coupled. The first was to 
obtain approximate Green function bounds taking advantage of the weakly coupled nature of the equations 
in the small-amplitude case to approximate by a superposition of solutions of scalar conservation laws; 
the second, by delicate pointwise interaction estimates between the approximate Green kernel and various 
algebraically decaying source terms, to close a nonlinear iteration and obtain the result. The "coupling" we 
mention refers to the fact that the Green function estimates blow up as amplitude goes to zero and the shock 
becomes more and more characteristic, whereas the source terms decay with amplitude; thus, the two effects 
must be delicately balanced to close the iteration and achieve a correct result. In the large-amplitude case, 
bounded away from the characteristic limit, the issues are somewhat different; namely, we do not have this 
"characteristic coupling" problem, but on the other hand we cannot as in |2II obtain approximate Green 
function bounds by asymptotic development in the amplitude. Likewise, for physical, partially parabolic 
viscosity, there are new difficulties associated with regularity and the need to gain derivatives. 

These new difficulties have largely been surmounted in the authors' previous work. In particular, (i) 
sharp Green function bounds have been obtained by Mascia-Zumbrun |22I in great generality using Laplace 
transform/stationary phase estimates, and (ii) global existence and sharp iJ* n LP estimates on the nonlinear 
residual v have been obtained by Raoofi [2Z| using the linearized estimates of [HI] , a key cancellation estimate 
of Liu, and a nonstandard energy estimate to control higher derivatives. These results are described in more 
detail in Sections|21and^ and we shall use them freely in our analysis. A major advantage of this "bootstrap" 
approach is that we need not close a nonlinear iteration, but only carry out a linear fixed-point argument 
to obtain our result. This was used in the previous work JH! to establish "nearly optimal" pointwise bounds 
in a relatively uncomplicated manner. As pointed out in |321Eli however, to obtain the full bounds of Liu 
requires estimates not only on v and as in but also characteristic derivatives {dt +a^dx)v as in j27| . 
To carry out these bounds requires an immense ammount of additional work, with consideration of numerous 
different cases, and accounts for most of the work of this paper. 

Remark 6. In the strictly parabolic case, the regularity requirement on the initial perturbation u = u — u 
may be significantly relaxed in Theorem^] from u S to u e (7*^+", a > 0, and the bound £^+0(l + |x|)^'^/^ 
imposed on |uo| alone, and not derivatives |9i^uo| and |9^ito|, by using parametrix bounds as described in 
|42[ 1401 1161 [T5] to establish short-time pointwise smoothing estimates. (That these bounds are necessary 
in the real viscosity case, on the other hand, may be clearly seen in the proof of Proposition [21) In the 
artificial viscosity case B = constant, which includes the case B = I treated by Liu the short-time 
existence/regularity theory becomes trivial j42l llfij . and we require no additional regularity beyond the 
single pointwise hypothesis |mo| < E + 0(1 -I- |a;|)~'^/^. 

Remark 7. Alternatively, we may subsume the strictly parabolic case (Bl) under the same analysis used 
to treat the partially symmetric case (A1)-(A3), substituting for parametrix estimates the more elementary 
Sobolev estimates of \6'6\ I37| . For, in this case, we may drop the requirement of symmetry of A±, which 
was used only to obtain a skew-symmetric K such that ^{KA -j- B)± > 0, whereas K = suffices when B 
is full rank. Since the An component is empty, the remaining structural conditions are trivially satisfied 
upon multiplying ^21 by a smooth symmetric positive definite A'^{u) such that ^A^B > (guaranteed by 
Lyapunov's Theorem). By a more careful accounting, taking advantage of associated simplifications, the 
results of HHIEIj with the exception of derivative bounds, may be obtained in the strictly parabohc case 
under regularity , matching the assumption of Goodman in the classic paper stability with respect to 
zero-mass perturbations, and the results of ^S] (for which are required bounds on dxu) under regularity 
H^. In the present case, for which pointwise derivative bounds are crucial to the argument, we require 

Plan of the paper. In Section[21 we recall the results of |33[I15| for our use. In Section^ we give the 
brief argument establishing Theorem ^ subject to certain integral estimates. In Section 0] we motivate the 
analysis to follow by reviewing a key cancellation estimate of Liu j25j upon which our analysis, and that of 
|27j . depends. In Sections El and Owe carry out the main work of the paper, establishing the deferred 
integral estimates. 
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2 Preliminaries 

We start by recalling some needed, known results. 

2.1 Profile estimates 

We first recall the profile analysis carried out in [331, generalizing results of "^T* in the strictly parabolic case. 
Profile u(x) satisfies the standing- wave ordinary differential equation (ODE) 

B{u)u' = F{u) ~ F{u_). (2.1) 

Considering the block structure of B, this can be written as: 

F\u^,u'')eeF\uL,uL') (2.2) 

and 

bi{u^y + 62(w")' = i^"(M^ m") - F^\uL,uL^). (2.3) 

Proposition 1. I]SS^J Given (H1)-(H3), 1^2. 3p determines a smooth r -dimensional manifold on which \2.!^) 
determines a nondegenerate ODE. Moreover, endstates u± are hyperbolic restpoints of this ODE, i.e., the 
coefficients of the linearized equations about u±, written in local coordinates, have no center subspace. In 
particular, under regularity (HO), 

DiDl{u\x) - u±) ^ 0{e~°'^''^) as X ^ ±oo, a > 0, < j < 10, i ^ 0,1. (2.4) 

2.2 Linearized equations and Green distribution bounds 

We next recall some linear theory from IJHIEIEI- Linearizing H1.2|l about u'^*(-), S■^, to be determined later, 
gives 

vt = Lv := -{Av)^ + [Bv^)^, (2.5) 

with 

B{x) := B(m** (x)), A{x)v := df{u^' {x))v - dB{u^' {x))vui' . (2.6) 
Denoting :— A(±oo), B^ :— B{±oo), and considering LemmaQ] it follows that 

|A(a;)-A-| = 0(6^"!^!), |B(a;) - B^l = 0(e-"l^l) (2.7) 

as X ^ — oo, for some positive rj. Similarly for and B^, as a; ^ +oo. Also \ A{x) - A^\ and \B{x) - B±| 
are bounded for all x. 

Define the (scalar) characteristic speeds af < ■ ■ ■ < (as above) to be the eigenvalues of A*, and the left 
and right (scalar) characteristic modes , r^ to be corresponding left and right eigenvectors, respectively 

(i.e., A^r^ = a^r^, etc.), normalized so that Ij^-r^ = Sj. and ■ — ^i- Following Kawashima |22| . 
define associated effective scalar diffusion rates (3f : j = 1, •■• ,n hy relation 



± n 



(at 



\ Pi 



diag L±B±i?±, (2.8) 



where := {if , . . . , l^Y, := {rf, . . . , r^) diagonalize . 
Assume for A and B the block structures: 



^^21 ^22 / ' \B2l B22 

Also, let a*{x), j — 1, . . . , (n — r) denote the eigenvalues of 

A^ := All — A12B22 B21, 
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with l*{x), r*{x) € R""'' associated left and right eigenvectors, normahzed so that Z**rj = S^. More generally, 
for an m*-fold eigenvalue, we choose (n — r) x m* blocks L* and R* of eigenvectors satisfying the dynamical 
normalization 



Lfd^R* = 0, 

along with the usual static normalization L*^Rj = S'^Im'] as shown in Lemma 4.9, this may always be 
achieved with bounded L*, R*. Associated with L*, R*, define extended, n x to* blocks 

c* ■= n* ( 5 

[oj' [-b-;b,,r; 

Eigenvalues a* and eigenmodes £*, TZ* correspond, respectively, to short-time hyperbolic characteristic 
speeds and modes of propagation for the reduced, hyperbolic part of degenerate system H1.2|l . 
Define local, mj x nij dissipation coefficients 

r]*ix) := -LfD,R*{x), j = 1, . . . , J < n - r, 

where 

D,{x) := [A21 - + ^,522^521 + B22d.,{B-^B2i) 

is an effective dissipation analogous to the effective diffusion predicted by formal, Chapman-Enskog expan- 
sion in the (dual) relaxation case. 

The Green distribution (fundamental solution) associated with (|2.5|l is defined by 

G{x,t-y):^e^'6y{x). (2.9) 

Recalling the standard notation errfn(z) := Jloo have the following pointwise description. 

Proposition 2. Under assumptions (Al)-(A3) [alt. (Bl)], (H0)-(H3), and (V), the Green distribution 
G{x,t\y) associated with the linearized evolution equations may be decomposed as 

G{x,t;y) = H + E + S + R, (2.10) 

where, for y < 0: 

J 

H{x,t;y) ^a*~l(a;)a*(y)7^*(a:)C;(2/,^)4^a.t(-y)/:f (y) 

7' (2-11) 
= ^7^*(x)0(e-''«*)5,_a.^(-y)/:f(y), 
i=i 

where the averaged convection rates a* = a*{x,t) in 1^2.11]) denote the time-averages over [0,i] of a*[x) along 
backward characteristic paths z* — z*{x,t) defined by 

dz* /dt = a*{z*), z*(t)=x, 
and the dissipation matrix — (*{x,t) G R'"^^'"^ is defined by the dissipative flow 

dC;/dt = -^*{z*)q, c;(o) = /™,; 

and 6x-a^t denotes Dirac distribution centered at x ~~ a"-t. 

Eix,t;y):=Y.%i^ ^^(2/'^)' (2-12) 

j — 1 |(5— (5* 



en 



{y,t) := I errfn I ^±S! - errfn | 1 IJ, (2.13) 



a. >0 
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Six,t;y) ■.= X{t>i} E r^lk\^^/3kt)-'/'e-^^-y-'^^'^'/'^^' 



a. <0 



>o 



>0, a - <0 



tfc >0, ti+>0 



and 



5± . , ^ a;± ^± , |y| / a; 



T/ie remainder R and its derivatives have the following hounds. 

n 

+ ^ O ((i + l)-i/2e-''-+ + e-"!"!) f-i/2e-(^-?/-«.-*)VMt 
fe=i 

+ E X{|.-,>|,|}0((i + 1)-!/ Vl/2),-(.-a-(t-|.K-|))VM*,-..+ 



>0, <0 

+ E X{|a^*|>|.|}0((t + l)-^/W^)e-(-4(*-l''K-|))VM*e- 



ar>0, a+>0 



t; t/) = E 0(e-''*)5,_a;t(-y) + 0(e-^(l--^l+*)) 
i=i 

n 

+ ^ O ((i + l)-i/2e-''^^ + e-"!"! + e""!^!) t-ie-(^-!/-«."*)VMt 
fe=i 

+ E X{|.-,|>H}0((i + l)-VVi)e-(-«7(*-l''K-|))VMte-'/^+ 

a^>0, o~<0 

+ E X{|„-,|>|,|}0((i + i)-i/vi)e-(-4(*-i.K-|))VM*,-..- 

a^>0, at>0 
.7 

y) = E 0(e-''*)(5,_a.t(-y) + 0(e-"(l--^l+*)) 

+ ^ O ((< + l)-ie-''^^ + e-"!"!) t'^t + i)i/2g-(x-!/-a,-t)VMt 
+ E X^|„-,|>|,|>0((i + l)-VV^)e-(— . 



k=l 

{t-\y/<^k\)flMt^-r,x+ 

a, >0,a, <0 



E X{|„-*|>|,|}0((t + i)-i/2,-i)e-(-a;(*-I.K-|))VMt^-,.- 



aj. >0,o+>0 
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Moreover, for \x — y\/t sufficiently large, \G\ < Ce~'''e~'^~^' Z*^*-' as in the strictly parabolic case. 

Setting G := S + B., so that G = H + E + G, we have the foUowing useful alternative bounds for G. 
Proposition 3 (|42|. | 33l Uti p. Under the assumptions of Proposition\^ G has the following bounds. 

J inax{0,|Q|-l} 

\d^,^Gix,t;y)\ < ^ 0(e-"*)a,^4-a*t(-y) + 0(e-''(l--^l+*)) 

J = l 13=0 



k=l 

E ,-1/2 -{x-aj{t-\y /a- \)f/Mt -nx + 

a^>0,a7<0 

I V V ^-1/2 -{x-a+{t-\y/a-\)f/Mt-rix-\ 



°-k >0,aj>0 



J 

\{dt + ar9.)G(a;, t-y)\<Y, 0{e-^')5^^-a-t{~y) + 0(e-"(l--«l+*)) 

j=i (2.20) 

(j^-:i/2 ^^-{x-y-a-tf /Mt^-rix+ _j_ ^ g-(a;-a- (t- ly/a" |)) VA«g-')2;+ ^ 



a. >0 



a->0 

+ C{t-^ + e-''l-l)(^i-l/2g-(x-y-a-t)VAftg-^.+ 



. >P V f-'i-/'2p-{x-a-(t-\y/a-\)f/Mt -rjx+ 

I V V ;f-l/2p-(a;-a+(t-|a/a-|))VAft 

+ 2^ ^{k,Tt|>|y|}^ e ^ ^ J' 

aj:>0, a+>0 

< |a| < 2, /or y < 0, and symmetrically for y > 0, for some rj, C, M > 0, where are as in Proposition 
> 0, denotes the positive/negative part of x, and indicator function X{\a^ t\>\y\} ^ \^k^\ — 1^1 
and otherwise. Moreover, all estimates are uniform in the supressed parameter 6^ . 

Remark 8. In the strictly parabolic case r — n, the hyperbolic part H is absent in the decomposition of 
G given in ifOH) . 

Remark 9. From (|2.13() . we obtain by straightforward calculation (see the bounds 

\e,iy,t)\<GY f-fn f^Verrfn 

\dte,{y,t)\<Gt-^/^ Y e-\y+-^*\"/'", 

a->Q 

\dyejiy,t)\<Ct-'^^ Y e-l^+'^'^*l'/^^* 

a->0 

\dyte,{y,t)\<Ct-' Y e-l^+^-^^lVM* 

a->0 

for y < 0, and symmetrically for y > 0. 



(2.21) 
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2.3 LP decay 

Finally, we recall the i/" n result of |37| established using the linearized bounds of and Hausdorff- 
Young's inequality together with an iJ* energy estimate and the cancellation estimate described in Section 
m below. 



Proposition 4 (| 37l [T^ ^. Under the conditions of Theorem^ there exists a unique, global solution u of 

^{5^)5, with |m± 



M.'4) . u e n . Moreover, for v := u - u^' - ip - ^{5*)S, with \mf\, \6^\ < CEq and 6{t) defined as 



m lO and J^TWjl . 

H-.t)\Lv <CEQ{\ + ty'^^^~p'>~-^, l<P<oo, (2.22) 

and 

H-,t)U.<CEo{l + t)-"^, (2.23) 

for some constant C > 0, 

3 Nonlinear analysis: Proof of Theorem [1] 

We now carry out the proof of Theorem ^ assuming certain integral estimates to be established in Sections 
El El andCl Let u be the solution of (fT^ guaranteed by Proposition H Following [221123 EZj, let |mi|, 
\5.^, \ < CEq he the imiquc solutions guaranteed by (I?ii) and the Implicit Function Theorem of 

/ + 00 
-oo _ I 

<0 aj>0 

or, equivalently, Jl^^ u{x,0) — J^^{u^' + (j)){x) dx. This determines the asymptotic state, by conservation 
of mass. 

Remark 10 (|2SIH2|)- In the case of Lax-type shock waves, u^{x) = u{x + 6), hence / dx = (u+ — u_) 

and (5* can be explicitly computed as the solution of a linear equation. 

Setting u{x,t) := u{x,t) — u^'{x), use Taylor's expansion around u^'{x) to find 

ut + iA{x)u)x - {B{x)ux)x = -(r(a;)(w, u))^ + Q{u, u^)^, (3.2) 
where r(a;)(u, w) = d'^f{u^*){u,u) — d'^B{u^*){u,u)ul* and 

Q{u,Ux) = 0(|u||u2;| + \u\^), 
with r± = r(±cx)). Define constant coefficients bf^ and F^j, to satisfy 

n n 

r±(r±,r±)=^r±,r±, B^rf^J2^rf. (3.3) 

4=1 1=1 

± _ ;,± „± .„ -p± 



Continuing, set v u — ip — ^^S{t), with S{t) to be defined later, and assuming S{0) — 0. Notice that. 



Then, of course, /3f — bf^ and := F 

ds 

by our choice of 6^ and diffusion waves (/jf 's, we have zero initial mass of v, i.e. 



v{x,0)dx = 0. (3.4) 

Replacing u with v + ip + ^^5{t) in H3.2|l ( computed at 5 — 5^,), and using the fact that satisfies 
the linear time independent equation Lv — 0, we obtain 

vt - Lv = $(x, t) + T{p,, V, —6{t)),— + 6{t), (3.5) 
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where 



(3.6) 



H^,v, ^S) = o(|.p + Mv\ + \v\\^6\ + M^6\ + \^s\' 

+ \iv + v + -Qg5)[v + V + + W + V + -Q^5\)- 

Furthermore, 

+ \iv + ^+^6U\iv" + ^ + ^5U (3.7) 

and ^{x,t) := -(ft - {A{x)ip)x + {B{x)ipa:)x - {r{x){(p,(p))a:. For $ we write 
^{x, i) = - [ipt + Aip.^ - Bip^x + T{ip, ip)x) 

= - E + iM^)v'rr), ^ iBix)^lr~), + (r(x)(^V-, ^V")), 

a~<0 

- E ^tr^ + (^(^)v'^t)^ - (BW.rth + (r(x)(^V+, <^V+)), (3-8) 

af>0 

Let us write a typical term of the first summation (a~ < 0) in the following form: 

V\ri + {A{x)ipWr)^ _ (B(x)^>r), + (r(x)(^V-, (^V")), 

= [(yl(a;) - A-)^Vr - {B{x) - B-)^lr- + (r(x) - r-)(^V-, ^V")]^ (3.9) 
+ ^Ir- + (iflA-r-) ~ is^JTr-) + ((^')2r-(rr, r,-)). 

Now we use the definition of Lp^ in (|1.10(l and the definition of coefficients hij and Fy-fc in H3.3|l to write the 
last part of (|3.9|) in the following form: 

^jr- + (<^;A-rr)-(^^^i3-r-) + {{^^-)lV- [r' ,rT)) 

Similar statements hold for > with minus signs replaced with plus signs. 
Applying Duhamel's principle, we obtain from H3.5|) 

/ + 00 
G{x,t\y)vQ{y)dy 
-OC 
ft p-i-OQ 



J J G{x,t- s;y)T{(p,v,—S)y{y,s)dyds, 

rt r+oo Q-5 

J J G{x,t-s:,y)<P{y,s)dyds + —S{t), 



(3.11) 
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where we have used the identity J^^ G{x, t\ y)^{y)dy = e^*^ = f|- and 5(0) = 0. Assuming 



/oo 
ei{y,t)vo{y)dy 
-OO 

ft r+oo Q-5 

/ e,{y,t- s)T{ip,v,—5)y{y,s)dyds (3.12) 

J-oo oo 

t /'+00 

/ ei{x,t- s;y)^{y,s)dyds, 

J-oo 



and using lprTT|l . (pDUl and G = H + E + Gwe obtain: 



/+00 
(H + G)ix,t;y)vo{y)dy 
-oo 

+ j {H + G){x,t-s;y)J^{^,v,—6)y{y,s)dyds (3.13) 

pt p+aa 

+ {H + G){x,t- s-y)<^>{y,s)dyds, 

Jq J -oo 

which we can clearly augment with derivative estimates through differentiation on both sides. In addition 



oo 



to v{x, t) and 5(t), we will keep track in our argument of Vx{x, i), {vt + a^Vx)^ and (j(t), the latter of which 
satisfies 

= - / dteiiy,t)vo{y)dy 

■/ — OO 

t r+oc Q-S 

/ dtei{y,t- s)T{ip,v,—S)y{y,s)dyds (3.14) 
J-oc 

t n+oo 

/ dteiiy,t- s)^{y,s)dyds, 

J-oo 

where we have taken advantage of the observation, apparent from (|2.21|) . that e{y, 0) = 0. 
Define 

C{t) := sup 



yfi<s<t -01(2/, s) + V'2(y, s) 

\vy{y,s)\ 



y,o<?<t s i/2(^^(y^s) +^2(2/, s)) +ip3{y,s) +ip4iy,s) 

\id, + afdy)v{y,s)\ (3.15) 



E 



sup . , _ . , 

a-<OM+>o^-°-"-' "^"'^^ + s)i/Vi' (y, s) + s-^/^{^i'^{y, s) + My, s)) + My, s) + My, 

\S{s)\ \S{s)\ 
sup \ ' + sup ' ^ ^' 



The desired estimates follow easily using the following integral estimates, to be established in Sections |S| 
and0 
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Lemma 1 (Linear estimates I.). If vo{y)dy = and |wo(y)| < Eo{l + \y\) Eq > 0, then 



/ + 00 
G{x,t;y)vo{y)dy < CEoil)i{x,t) 
-OO 

+ 00 

d,G{x,t;y)voiy)dy < CEoir^/^ + e-^^^^)^^,*) 

-OO 

idt + afd,)G{x,t;y)vo{y)dy < CEo(t-^i4'^{x,t) +t-'/m'^{x,t) + e-^\'^\Mx,t)) (3.16) 

-C!0 

/+00 
ei{y,t)vo{y)dy < CEo{l + t)-'/^ 
-OO 

/ + 00 
dtei{y,t)vo{y)dy < CEo{l + t)-^/^. 
-co 

Lemma 2 (Linear estimates II.). If\vQ{x)\, \dxVQ{x)\, \d^VQ{x)\ < EQ{l + \x\)~i , Eq > 0, then, for some 
0>O, 

/+00 
H{x,t;y)vo{y)dy < CEoe-^\l + \x\)-^ < CEoMx,t), 
-OO 

/ + OC' 
H,{x,t;y)vo{y)dy<CEoe-'^\l + \x\)-i <CEo{l+t)-'/^Mx,t), (3.17) 
-co 

/+00 
{dt + afd^)H^{x,t;y)vo{y)dy < CE^e-^^l + \x\)—- < CEo{l + ty^/^M^,*)- 
-OO 

Lemma 3 (Nonlinear estimates L). If \\v{y,s)\\Loo < CEo{l + s)'^/^, \\v{y,s)\\m < CEo{l + sy^/^ , 
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and (^(t) < +00 fsee \,y.l5\) . then 

rt /•+00 



10 J~oo 



G{x,t - s;y)T{ip,v, —5)ydyds < CEQ({t)[ijji{x,t) +^2(2;,*)] 
G{x, t-s;y)^y, s)dyds < CEoipi {x, t) 



'0 J-00 



Gx{x, t - s; y)J^{(p, v, -^S)ydyds 



'as 



{dt + a±9.) 



/O J-00 
ft p+oa 



10 J-00 

< CEoCit) [t-^^\Mx, t) + M^, t)) + Mx, t) + t)] 

ft p+oa 

G^{x,t- s-y)<^{y,s)dyds < CEot-^/^Mx^t) 

G(x, t- s; y)T{ip, v, —S)ydyds 

< CEoCit) [t-\l + t)'^^iPi^^ix, t) + t-'/^{^i^^{x, t) + ^2{x.. t)) + ^3(2^, t) + 4>i{x, t)] 

t p-\-OQ 

G{x,t~ s-,y)^{y, s)dyds 



J -00 
< CEo 

t /'+00 



'0 J~oo 



'0 J-00 



dye,{y, t - s)T{^, V, —6)ydyds < GEoC{t){l + t)-""^ 
dyte^{y,t - s)T{ip,v, —S)ydyds < CEoC{t)il + t)-^ 



'0 J-00 



e,(y, t - s)$(2/, s)dyds < GEoil + t)-^'^ 
dte,{y,t~ s)<^{y,s)dyds < CEa{l + t)-\ 



J-00 



(3.18) 



Lemma 4 (Nonlinear estimates II.). U \\viy, s)\\l^ < CEo{l + s)-^/^ , \\v{y, s)\\Hr. < GEoil + s)-'^/'^ , 
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and (^(t) < +00, then 



H{x,t-s- y)T{ip, V, -^5)ydyds < CEoCit) [V'l {x,t) + ^2 {x, t)] 



/ / H{x,t- s;y)^{y,s)dyds 

Jo J -00 



J Hx{x,t- s]y)T{ip,v,—5)ydyds 



{dt + afd^) 



{dt + a±9.) 



< CEoCit) \t-^'^{Mx, t) + t)) + ^3(2;, t) + iji{x, t)] 

t p+oa 

H,{x,t-s;y)^y,s)dyds < CE^^t-^'^i^Ji ^ i,2){x,t) 



t /' + 00 



H{x, t ~ s; y)T{ip, v, —6)ydyds 

< CEoCit) [t-\l + t)'^^iP'{^ix, t) + t-'/^{^i^^{x, t) + Mx, t)) + Mx, t) + 4>i{x, t)] 
/ / H{x,t~ s;y)^{y,s)dyds 

Jo J -oc 

< CEo \t-\l + ty/*^Pi'^{x, t) + t-^^\'4^i'^{x, t) + ij2{x, t)) 



(3.19) 



Proof of Theorem U We prove Theorem ^ directly from the integral representations (|3.12(l and l|3.13|l . 

augmented with similar representations for 5i{t), Vx{x,t), and {dt + a^dx)v{x,t), obtained through direct 
differentiation of H3.12|l and (|3.13(l . Recalling the definition of our iteration variable (^{t) in H3.15|l . our goal 
will be to employ the estimates of Lemmas ^ 01 to establish the inequality 



c{t)<cEo{i+at)), 



(3.20) 



where Eq is precisely as in Theorem ^ Choosing, then, Eq < and noting that the possibility of ({t) 
jumping discontinuously from a finite value to an infinite value is precluded by short-time theory (see Remark 
3.6 in n^)-'^, we will be able to conclude 



CEn 



1-CEo 



< 1. 



(3.21) 



The estimates of Theorem ^ follow immediately from (|3.21|l and H3.15|) . 

We begin with a careful consideration of the integral represention for v{t, x), (|3.13|) . which we will separate 
for clarity into terms that arise in the case of strict parabolicity — involving G — and the additional terms 
arising from the relaxation of strict parabolicity — involving H. For the terms involving G, we estimate 



+00 



G{x,t;y)vo{y)dy 



+00 



G{x, t - s; y)T{ip, v, -^5)yiy, s)dyds 



~d5 



+00 



G{x,t- s;y)^{y,s)dyds 

10 J -00 

< GEo^Ji {x,t) + CEoCit) [^1 {x,t)+Mx,t)] + CEo^i (x, t) , 

where the estimates follow respectively from the first estimate of Lemma ^ the first estimate of Lemma 01 
and the second estimate of Lemma [31 (though the the first and last estimates are the same, we write both 

^Note also that the pointwise bounds of jlSl directly imply that is finite for each given t. 
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for clarity). For the terms in (|3.13|) involving i7, we have similarly 



+ C30 



H{x,t; y)vo{y)dy 



t /' + 00 



H{x, t - s- y)T{(p, V, -^S)y{y, s)dyds 



^-oo 



/ / H{x,t- s;y)^{y,s)dyds 

Jo J-oo 

< CEM^, t) + CEoat)[Mx, t) + M^, t)] + CEoMx, t), 



for which the estimates follow respectively from the first estimate of Lemma |21 the first estimate of Lemma 
21 and the second estimate of Lemma 0] Combining these estimates, we conclude 



\v{x, t)\ < CEoiM^, t) + t)] + CEoCit)[Mx, t) + t% 
which can be rearranged as 



li)l{x,t) +1p2{x,t) 

Keeping in mind that is a nondecreasing function of i, we have at last 

Hv^s)\ 



sup - , > , , - 

vfi<s<t vi{y,s) + 'ip2{y,s) 



< cEo{i + at))- 



(3.22) 



Proceeding similarly for each of the expressions 6i{t), 6i{t), Vx{x,t), and {dt + dx)v{x,t), we can bound 
each summand in the definition of C,{t) in precisely the same way by CEq{\ + C,{t)). We conclude the sought 
estimate 

C(t) <C£;o(i + CW)- 

As discussed in H3.2()|l and (|3.21(l . we conclude C,{t) < 1, and from this the estimates of Theorem ^ □ 



4 Liu's cancellation estimate 

Before carrying out the deferred integral estimates, we revisit a key estimate of Liu that at once deter- 
mines the ultimate rate of decay and motivates the analysis to follow. Consider the illustrative convolution 



u{x, t) 



+00 



+00 



■j{x -y.t- s){K(y ~ s, sf)y dy ds 



,{x-y,t- s)K{y - s, sf dy ds, 



(4.1) 



g{x,t) = K{x,t) = (47ri)^2eT^, similar to quadratic interaction integrals arising through the integration 
of scattering terms against diffusion waves (see Remark II l|l . If we replace the integrands in l|4.1|l by their 
absolute value, we obtain (see ^H]) the sharp estimate 



\u{x,t)\<C{g{xAt)+9{x-tAt)) +Cx{^<,<t_^}(a; ^{t-x)-^). 



(4.2) 



By taking account of cancellation, however, we may obtain the following stronger bound (also sharp) 
pointed out by Liu We follow the notation of Raoofi and also the proof, based on integration by 
parts in the characteristic direction, which abstracts from the more concrete calculation of Liu the central 
ideas that will be of use here. 

Proposition 5 ( [251 137] \ For u{x,t) defined in j/^. j| ) and t>l, 

\u{x, t)\ < Ct--^ (gix, 8t) + gix ~ t, 8t)) + Cx^^t<x<t^^t] {t-Ht - x)-' + x'^) , (4.3) 

where x stands for the indicator function, and C is a constant independent of t and x. 
The same result holds if in is replaced with Kx ■ 
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This was first pointed out in 

Proof of Proposition\^ As ^ K^, the proof will be stated only for . It is straightforward to verify 
that the same argument works for at every step. We first state a simple lemma. 

Lemma 5. If < s < ^/t, then e J* < Ce~s^ with C independent oft,s and x. 
Proof. The statement of the lemma is equivalent to 

-{X±sf -2 



At 



< 



8t 



D 



for some D, which (after some calculation) in its turn is equivalent to 



{x ± 2s)2 - 2s2 > -8Dt, which holds for D > i, since < t. 

The argument relies on the following simple properties of the heat kernel g. 

9(x - y, t)g{y, t')dy = g{x, t + t') 

\gx{x,t)\ <Ct-^g{x,2t) 
\gt{x,t)\ < Ct-^g{x,2t) 

\g{x,t)\<Ct-^. 



Rewriting (|4.1|l . we have 



u{x, t) 



t p+oo 



g[x -y,t- s) {g{y - s, sf)y dy ds 



J-oo 



g{x -y,t- s){g{y - s, sf)y dy ds 



Vt 
t 



g{x -y,t~ s){g{y - s, s)^)^ dy ds 



OO 
+ CX2 



' t — y/t J — OO 

=: I + II + III. 



g{x -y,t- s){g{y - s, sf)y dy ds 



(/) and (///) are easy to estimate: 



g{x -y,t~ s){g{y - s, sf)y dy ds 



J -OO 







gy{x -y,t~ s)g{y - s, sf dy ds 



By ()4.5|l and (|4.7() . the above is less than or equal to 



c 



[t-s] --S ^g{x- y,2{t- s))g{y - s,2s)dyds 



which, by (|4.4(l . is less than or equal to C J^{t — s) ■^g{x — s, 2t) ds. Now, using LemmalHl the above 
is 



□ 



(4.4) 

(4.5) 
(4.6) 
(4.7) 



(4.8) 



<Cg{x,At) I (t- sy^s^^ ds < Ct'^^gix^At) I s^"^ ds 

< Ct-ig{x,At). 



(4.9) 
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Part (///) in (|4.8|l can be handled similarly. 
The more difficult part is part (//) of (|4.8I) : 

//=/ / g{x-y,t~s){g{y~s,s)'^)ydyds. 

JVi J -oo 

In order to estimate (//), let us write g — g(x, r). Then we have 

g{x -y,t- s){g{y ~ s, sf)y = {g{x -y,t- s)g{y - s, sf)s 

- 9t{x - y,t - s)g{y - s, sf 
+ 9{x-y,t- s){g'^)r{y'- s,s). 

We will do the estimates piece by piece. The first part of H4.11|l can be decomposed as 

/•t-Vt f+C 



(4.10) 



(4.11) 



{gix -y,t- s)g{y - s, s)2)s dy ds 



+ 00 



g{x - y,Vi) g{y - t + Vi,t - -Jtf dy 



+ 00 



g{x -y,t~Vi) g(y - Vt, Vif dy. 



(4.12) 
(4.13) 
(4.14) 



Using 14. 4() and 14. 7|) . it follows that 



+00 



g{x -y,Vt)g{y -t + Vi,t- Vtf dy < Ct ig{x-t + Vi,t), 



and 



+00 



g{x -y,t~ Vt) g{y - \/t, Vtf dy < Ct « g{x - Vt,t), 



but, by Lemma|Sl g{x — \/t, t) < g{x, 2t) and g{x ~ t + \/t, t) < g{x ~ t, 2t). These terms fit in the right hand 

side of E2J). 

For the other parts in H4.11|l . we use (|4.6() to obtain 

rt-Vi n+o 



/ / 9t{x- y,t- s)g{y ~ s,sf dyds 
J \ft J —00 

s~^{t - s)'^g{x - s, 2t)ds 



Vt 



< 



s 2(t-s) ^g{x-s,2t)ds 

+ S 2(t_s) ^g{x - s,2t)ds 

=:A + B. 



(4.15) 



If X < then 

A<Ct-^g{x~ Vi,2t) I s-ids 

< Cr^ g{x ~ ^/^, 2t) 

< Cr^g{x,At). 
Similarly when x > t — y/t, we obtain A < Ct^^g{x — t, At). 



(4.16) 
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Now ior ^/t < X < t — ^/i we have 



s ^{t-s) ^g{x - s,2t)ds 

t_ 

< Cr^ / ' s^^g{x - s,2t)ds 

r— r- 

= Cr^ ' s^ig{x- s,2t)ds + r^ j s^i g{x - s,2t)ds 
JVt J§ 

X n 1 1 1 /"^ 

<Ct"\g(-.2i) / s'^ds + r^x~^ I g{x - s,2t)ds 

2 Jvt J% 



(4.17) 



< C {t-^g{xAt) +t-'^x''^^ 



also acceptable, as t ^x < x 2 for \ft < x <t — \fi. 

Part B in (|4.15l) can be estimated similarly. We carry it out briefly only for | < a; < i — \ft. Let 
^ --t- Vt- x. Then 

B<Cr^\ I +J j {t- s)-^g{x- s,2t)ds 



< Ct-Ht x)-' + Ct-ig{ * ^ ^ 2t), 



to which we apply Lemma |S1 

There remains the last part of 14.11II . i.e., 



t-Vt r+ 



l{x -y,t- s){g^)riy - s, s) dyds 



which can easily be shown to be less than or equal to 



C s 2 g(^x - s, 2t) ds. 



li X < y/t, then (|4.20() is of the order 

Cg{x - Vi,2t) / s^i ds 

< Cr^ g{x- Vt, 2t) 
< Cri gix,At). 

For a; > t ~ y/t we use a similar method. 

For Vi < x < t - y/t, we use a similar method to those used for the previous cases to get 

(|On|l < c(rig{x,4t)+x- 



(4.18) 



(4.19) 



(4.20) 



This completes the proof. 



□ 
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Remark 11. Replacing K{y~at, t) by K{y, t) in H4.1|l . and g{x — y, t) by G{x^ t; y) we obtain bounds similar 
to H4.3|l (with appropriate modifications due to different speeds) provided k and G satisfy 

\g{x,t;y)\<Cg{x-y-at,pt), (4.21) 

\gy{x,t-,y)\ < Ct-^g{x-y-at,2(3t), (4.22) 

\gtix,t;y)\ < Ct-'gix ~ y - at,2pt), (4.23) 

\K{x,t)\ <Cg{x-bt,pt), (4.24) 

\Ky{x,t)\ <Ct-^g{x-bt,2l3t), (4.25) 
\Kt{x, t- y)\ < Cr^g{x - bt, 2(3t), (4.26) 

for some a ^ b^ and some constants C, /3 > 0, as hold in particular for k a single diffusion wave and Q a single 
component of the scattering term S in the decomposition G = E + H + S + R. Comparing with (|1.14l) . we 
see that these are the rate-determing terms. 

More generally |42j . we may expect analogous cancellation whenever k and G have the property that 
they decay faster along characteristic directions dg + a^dy [resp. dt + a^dx] than dy or ds [resp. or dt]. 
This is the underlying principle guiding our analysis, and that of [27) . 

Remark 12. Reviewing (|4.3|l more carefully, we see that the signal u is slightly more concentrated near 
the characteristic direction dx/dt = of the propagator g than in the direction dx/dt = 1 of the source 
Ky, decaying as {x + t^)^^ rather than t~^{\t — x\ + t4)~i. Systematic, characteristic-by-characteristic 
bookkeeping taking account of this difference leads to the slightly refined bounds of |27]; these sum, of 
course, to the simpler modulus bounds presented here. 

5 Linear integral Estimates 

It remains to establish the deferred integral estimates used in Sectional We begin, in this section, with the 
linear integral estimates of Lemmas ^ and El 

Proof of Lemma ^ The first, fourth, and fifth estimates of Lemma ^ have been established in [Hj. The 
second and third follow in a similar fashion from the estimates of Proposition |3| □ 
Proof of Lemma^ Looking at H2.11|l . we notice that in order to estimate H{x, t; y)vo{y)dy it suffices 



f -t-oo 

-oo 

to estimate 



TZ*{x)0{e"'^°')S.,^a>t{-y)C*\y)vo{y)dy < CEoe-'^«'vo{d*t - x) 



Here we used the crude inequality 



l + \a + b\ - 1 



<CE„e-^''\l + \d*t-x\)-'^ (5.1) 

< CEoe-"»'il + \x\y^l + \d*t\)^ 

< CEoe-^{l + \x\rl. 



and the fact that a* TZ* and £** are bounded. Observing that e < C{N){1 + t) for any N, and 



(1 + |x — at\) < C(l + t){l + \x\), we can bound the righthand side of H5.1|l in turn by CEqc ™ ^'1(2;,^), 
giving H3.17|l fi^. Estimates 1)3. 17|l fill and H3.17|l fiii^ are obtained similarly. □ 

□ 

6 Nonlinear integral Estimates I 

In this section, we carry out the main work of the paper, establishing the nonlinear integral estimates of 
Proposition El 
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Proof of Lemma 121 Under the assumption that C(^) is bounded, we have the estimates 



\vix,t)\<at) 



(x,t), 



\d.,v{x, t)\ < C{t) [t-^/^{^l + ij2) +^3+ ^-4] {X, t) 

\{dt + ap.,)v{x, t)\ < at) [t-\l + + + ^2) + V-s + ^4] (x, t). 

In the analysis that foUows, we will omit C,{t) from our estimates in most cases and focus only on the terms 
with a given form, the template. 

We observe at the outset that in proving estimates of form ipiix^t), we will frequently make use of the 
inequality 



(1 + 



-3/4, 



^ <C(l + |x-a±t|+t-i/2)-3/2. 



_ . , . (6.1) 

In the case \x — a/^t\ < KVt, for some constant K, this inequality is immediate. On the other hand, for 
\x — o-kA — KVt, we observe 



(1 + i)-3/4e-^^^^ ^ (1 + t)-y^\x - a±ir3/2i3/4 



\x - a±i|3/2 



^3/4 



< Cl(l + t)-3/4i3/4|^ _ ^±^|_3/2g-^i^J^ < _^ 1^ „ ^±^1 ^ il/2)-3/2^ 

where the seeming blow-up as |x — a^t| — > is controlled by the size of y/t, which must be smaller than 
Ix-aftl. 

Proof of ()3.18f i')^. For the first estimate in Lemma El we begin by considering the nonlinearities 
s)'^)y and (v(y, s)vy{y, s))y. Here and in the remaining cases, the analyses of the convection, reflection, 
and transmission contributions to the Green's kernel dyG{x, t; y) are similar, and we provide full details only 
for the case of convection. We observe that the contribution 

J max{0,|Q|-l} 

^ J2 0{e-^')dl5,.-,n{-y) 
3=1 P=o 

is similar to, though less singular than, terms arising in dyH{x,t;y) (see 12.1 J|) . and can be analyzed as in 
the proof of Lemma 01 Finally, we remark that the contribution 

has no effect on the iteration. 

^S.lt^ i)). term one. We first consider integration of our convecting Green's kernel against the nonlinearity 
s-i/2(^]^ + s)^'^/'^ilji{y, s). In this case, we have integrals of the form 



[t - s)-^e- ^ ° (1 + \y- als\ + s^'^Y"^'^ s-^'^{l + s)-^'^dyds, 



t - s)"ie" 

/O J-00 

with a similar integral for y > 0. Proceeding as in |15) . we write 

X ~y-aj{t- s) = {x- aj{t - s) - a^s) - (y- a^s), 
from which we have the estimate 



(6.2) 



(6.3) 



(1 + \y~ a-s\ + si/2)-7 

(x-y-a-(t-s))2 (ic - a " ( t - a ) - a " a) ^ 

e"' ^TTrrT) e (1 + |y - a^s| + s^'^)"'^ 

(x-y-aJ(t-3))2 

_^ g Mjt^ (1 + |y - a^TsI + \x-aj(t- s)- a'^s\ + s^/^)""*" 



< C 



(6.4) 
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where 7 = 3/2. (Here and in future cases, we will state estimates in terms of a parameter 7 for general 
applicability.) For the first estimate in H6.4|) . we have integrals 

(t - 3)-^-'^ \^,it-s) " e~" '*^('--) (1 + |y - s| + si/2)-3/2^-i/2^^ ^ sy^/^dyds. (6.5) 



^-00 



We have three cases to consider: < < a~, < a~ < 0, and a~ < < 0. We will proceed by 
analyzing the second of these in detail and observing that no qualitatively new calculations arise in the 
remaining two. For this second case and for > IflfTlt, we write 



X - {t- s)- s^{x- t) - (a^- - )(i - s). 



(6.6) 



for which we observe that there is no cancellation between summands. We immediately obtain an estimate 
on H6.5|) by 

^ ^ r -f I / 



Cit-'e 



(l + sl/2)-l/2s-l/2(l + s)-l/4^5 



(6-7) 

t/2 

< C(l + t)~Mn(e + t)e~ '"'^'" , 
which is sufficient by H6.1|l . For |a;| < |a~|t, we write 

X- a'j{t- s) - a^s = {x - ajt) - {a'^ - aj)s, (6.8) 

for which we observe that there is no cancellation between summands, and proceeding as in (|6.7|l . we obtain 
an estimate by 

C(l + t)-^ ln(e + t)e- '" ^'^ " , 

which again is sufficient. We remark that this concludes the analysis in the case aj = ajT, so from here on 
we may take ojT < aJ and the case \o.J\t < \x\ < |a^|i. In this case, and for s € [0,t/2], we observe through 
(j6.8(l the inequality 



(x-a-(t-s)-a-a) 



' (1 + \y- a-s\ + si/2)-3/2(i + 5)-T 



< C 



(x-a t) 



e ^{1 + \y~a~s\ + s^/^)-^/^{l + s)-'' 



(6.9) 



(x-a7(t-s)-a. s) 



+ e mT^) {l + \y- a-s\ + |x - a7i|i/2)-3/2(^ ^ |^ _ aj^l)"'^ 



with 7 = 3/4. For the first estimate in H6.9|l . we proceed similarly as in H6.7|l . while for the second, we have, 
upon integration of the e-kernel in y, an estimate by 



Cit-V2(i + 1^ _ „^-^|i/2)-3/2(^ ^ _ ajt\)-^/^ 
<Cil + \x-aJt\r^/^, 



t/2 



(x-a, (t-s)- 



M(t-.) 



where we have used in this last inequality that aJ =/= ■ We observe that this is clearly sufficient in the case 
|a; — ajt\ > ^/t, whereas in the case \x — ajt\ < ^/i, decay of kernel type exp((a:; — ajt)'^ /{Lt j} is immediate, 
for L sufficiently large, and so we only require decay at rate which is straightforward. For s £ [t/2, t], 

we observe through (|6.6I) the inequality 



{t - sy^e mj^) 



< C 



(6.10) 
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where 7=1. For the second estimate in (|6.1UI) . we proceed similarly as in (|6.7|) . while for the first we have, 
upon integration in y of the nonlinearity, an estimate by 



< C{l + t)-^'^\x-a-t\-^. 



t {^_„-(t-»)-a~s)2 

e 

t/2 



M(t-s) 



ds 



Recalling that we are currently working in the case t > \x\/\af. |, we see that this final estimate is sufficient 
for the case \x — a^t\ > \ft. For the second term in (|6.4|l . we have integrals of the form 



, „ 

/O -/-oo 



(x-B-a. (t-s)) 



(t - s)-^e s^rfc^) (1 + |j/ - a^s| + |x - a7(t - s) - a^Ts] + s^l'^Y^I'^ s-^l'^{\ + sY^l^ds. (6.11) 



We continue to focus on the case a^. < a - . For the cases \x\ > |aj, \t and |a;| < \aj \t, respectively, we have 
no cancellation between the summands in (|6.t)|) and (|6.8|) and consequently we obtain estimates 



C(l + (l + \x- a-tiy-'r' + il + \x- a;t\) 



3/2 



-3/2 



(6.12) 



For the case \aj \t < \x\ < ja^, \t, we divide the analysis into sub-intervals s E [0,t/2] and s E [t/2,t]. For 
s E [0,t/2], we observe through 16.8|l the inequality 



/ \ -3/2 

[l + \y- flfe s| + \x^ aj{t - s) - a^s\ + s^/^j (l + g) 



< C 



(l + |a;-a7t|+.si/2)-3/2(i + 5)-7 



(6.13) 



+ {l + \y- a-s\ + \x- a-{t - s) - a-s\ + \x - aVt|i/2)-3/2(^ _^ |^ _ ^-^^y 



for 7 = 3/4. For the first estimate in H6.13|l . we proceed as in H6.12|) . while for the second we have, upon 
integration of the kernel, an estimate on H6.11|l by 

/•t/2 

Cit^i/2(1 + \x- a-t|)-3/4 / ^^-i/2(^ ^ ^)i/2(^ ^ 1^ _ -(^ _ ^) _ ^-^1 ^ 1^ _ aTt|i/2)-3/2ds 

Jo 

<C{l+t)-^'^{l + \x-aJt\)-\ 

For s E [t/2, i], we observe through (|6.6|) the inequality 

{t - sy\l + \x~ a-{t -s)~ a^s\ + s^l^)-^'^ 

<c\\x- alt\-\\ + \x~ aj{t - s) - a' s\ + ,sV2)-3/2 ^^^^^^ 



+ it- s)-\l + \x -aj{t- s) - a^s\ + \x - ajt\ + s^/^) 



-3/2 



For the second estimate in H6.14|l . we proceed similarly as in ()6.12|) . while for the first we have, upon 
integration of the kernel, an estimate by 

C2\x ~ a^t\-H-^^\l + t)-i/4 /* {t - s)i/2(i + \x-a;{t-s)^ a^s\ + t^'^y^'^ds 

Jt/2 

< C\x - a-t\-\l + t)-i/4(i + 

which is sufficient for t > |a::|/|a~|. 

yj.lt^ i)). term two. We next consider integration against the nonlinearity s~^/'^{\ + s)~^/^-i/'2(j/, t), for 
which we have integrals of the form 

f f {t - s)-ie-''"''-°*-^r" (1 + + s)-''\l + |y| + s)-^'\l + \y\ + s^'^Y^/^dyds, (6.15) 

^0 J~\a-^ \s 
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wherein we have observed that for y G l^i |s,0], s decay yields also decay in y. We first observe an 
immediate time decay estimate by 

C,r^ / (f + + sl/2)-l/2^^ ^ C2t-^l\\ + t)-3/4(l + <l/2)-l/2 / „ sY^l^ds 

Jo Jt/2 (6.f6) 

< C{l + t)-Hn{e + t). 

In order to determine estimates in space as well, we observe the inequality 

(x-»-a7(t-s))2 



(1 + \y\r'^'{\y\ + s)-'/Hi + \y\ + + + sV2)-i/2 



< c 



e 



= °" e-'^ -(-^)" (I + + s)-i/2(I + |y| + sy'/^l + \y\ + s^'^T^I^ 

+ e-'"' "mc^-s)'" (1 + |a; - a-(t - - a-(t - + s)-^'^ 

X (f + |x - a-{t - s)\ + s)-3/4(i + _ _ s)| + si/2)-i/2 

For the first estimate in Ht).I7|l . we have integrals 

/ / it - s)~^e " e 

Jo i-lark (6-18) 
X (I + \y\)-'^'i\y\ + + \y\ + 5)-3/4(l + |y| + si/2)-i/2rfyd,. 

We have two cases to consider here, aJ<Q and aj > 0, of which we focus on the former. In this case, for 
> 1^7 1*' there is no cancellation between x — ajt and ajs, and so proceeding as in H6.16|l . we obtain an 
estimate by 



C(l + t)-^ ln(e + t)e-^—^^. 

For \x\ < \a~\t, we divide the analysis into two cases, s £ [0,i/2] and s E [t/2,t]. For s E [0,t/2], we observe 
the inequality 

e"'^ M^f Qy\ + + + s)-3/4(i + Ij^I + 5l/2)-l/2 

< c 



e ^ (|y| + + |y| + 5)-3/4(i + |y| + ^,i/2)-i/2 (6.19) 

+ e"'^ (s + |x - a7i|)"'/'(l + s + |a; - a7t|)"'/^(l + s'/' + |a: - a7t|i/2)-i/2 

For the first estimate in H6.19|l . we proceed similarly as in H6.16|l . while for the second we have, upon 
integration of the e kernel, an estimate by 

ft/2 (x-„-(t-s))2 

Cit"i/2(|a; - aTt\)-^/^{l + \x - ajtl)-^ / e ds 



< 



C{\x-aJt\)-^^\l + \x-a-t\)-\ 



where the seeming blow-up as |a; — <| can be eliminated by proceeding alternatively for |a; — t\ 
bounded. For s E [t/2,t], we have 

(x-„-(t-s))2 (x-a-(t-s))2 

(t-s)~^e ^^^<^> < C|a;|~^e 2M(t-s) ^ (6.20) 
7 = 1/2, for which we have an estimate by 

C2|xrl/2i-l/2(i^^)-3/4(^^^l/2)-l/2 e-'-^^^^T^ds<C\x\-'/^il+tr\ 

Jt/2 
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where the seeming blow-up as can be ehminated by an alternative calculation in the case of |a;| 

bounded. In the current case of \x\ < \aj\t, this final estimate is bounded by ip2{x,t). For the second 
estimate in (|6.17|l . we have integrals 



{t - sy'e ^TTT^ (1 + \x-a-{t~ s)\y 



-1/2 



'0 "'-la-Is " (6-21) 

X {\x - aj{t - s)\ + + \x- a-{t - s)\ + s)-'"\l + \x- aj{t - s)\ + s^'^y^'^dyds. 

Focusing again on the case a~ < 0, we observe that for |x| > |aj|i, there is no cancellation between x — ajt 
and a"s, and consequently that we have an estimate, upon integration of the kernel, by 



C,t''/^il + \x~aTt\)-y' 



t/2 



a-{t~s)\ + s)-^^^ds 



+ C2{l + \x^ a^t\)-''H-^'^ / [t - sr^'^ds 

Jt/2 

<Cil + \x-aTt\)-y\ 



(6.22) 



which, for \x ~ a t\ > ^/t decays faster than the claimed estimates. For the case \x — a ■ we require 

only t^^^'^ decay, which is clear from (|6.1t)|l . For < \aj\t, we divide the analysis into cases, s £ [0,t/2] 
and s e [t/2,t]. For s € [0, i/2], we observe the inequality 

(1 + |x - ajit - s)\)-'/'{\x a-{t s)\ + s)-^/^ 
X (1 + |x - a]{t ~ s)\ + s)-3/4(i + _ aT{t - s)\ + si/2)-i/2 



< C 



{l + \x- a-t\)-^/^{\x - a-t\ + .s)-i/2(l + \x - ay\ + s)-^'^{l + \x- a7<| + s^''')-^''' 



(6.23) 



+ (1 + Ix - a-{t - s)r i/2(|x - a]t\ + s)-^'^{l + \x- ajtl + sY^'^l + ajt\^^^ + s^'^)-^''^ 
For the first estimate in H6.23|l . we obtain an estimate, upon integration of the kernel, by 

Cit-^'\l + \x- a-t\)-^ / {\x - a]t\ + s)-^'^{l + |x - aTt\ + sY^'Hs < Cr^'^il + \x- aJt\)-^'\ 

Jo 

while for the second we obtain an estimate by 

.t/2 

Ciri/2(1 + \x~ a7t|)~3/2 / -aJit- s)|)-i/2ds < C(l + \x ~ ajtiy^'/^ 



For s E [t/2, i], we write 

(i_,)-l/2(l+|^_«-(^_,)|)-l/2 



< c 



\x\-'^^l + \x~aT{t~ s)|)-i/2 + _ s)-i/2(l + \x-aj{t^ .s)\ + \x\)-'/^ 



(6.24) 



For the first estimate in H6.24|l . we obtain an estimate by 

C2|x|-l/2t-l/2(l + t)-3/4(i ^ il/2)-l/2 (1 + 1^ _ aT(t - s)\)-'/^ds < C(l + t)-'\xr'/^ 

Jt/2 

which for |a::| bounded away from is bounded by ^2{x, t), while for the second estimate in (|6.14() . we obtain 
an estimate by 

C2(l + |x|)-l/2t-l/2(l + i)-3/4(i + il/2)-l/2 f' ,)-l/2^,. < ^ {I + \x\)-^'\l + t)-\ 

Jt/2 
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yj.ll^ i)). term three. We next consider integration against the nonlinearity (1 + s) ''l^' (which arises, 



for example, from the term (^t-(5)^), for which, in the case of our convection Green's kernel estimate, we 



have integrals of the form 



(x-y-a (t-s))2 

{t-s)-^e "(t^^o {1 + sY^er'^^y^dyds. 



(6.25) 



/O J -ao 

First, we observe a straightforward time decay estimate by 



t/2 



t-1 



t/2 



[t- s)-^ds + C3{l + t)-^ / {t-s)-^/^ds 



Cit'^ / {l + s)-^ds + C2{l + t)-^ 
Ja 

< C{l + t)-'^\n{e + t). 
In order to determine estimates in space as well, we observe the inequality 



(6.26) 



(x-»-a. {t^3)r 



< 



v\y\ 



(x-a7(t-s))2 



(x-j,-a (t-a))-' 

M(t-s) „-')l|^-Qj (*-s)L-')2|y| 



(6.27) 



For the first estimate in H6.27|l . we have integrals 

rt rO (x-„-(t-.))2 



{t-sy^e (1 + s)-ie-''l?'ldyds 



(6.28) 



We have two cases to consider, < and aj > 0, of which we focus on the former. For \x\ > \aj \t, 
there is no cancellation between x — ajt and aJ s, and proceeding almost precisely as in (|6.26|l . we obtain 
an estimate by 

C(l + t)-^ ln(e + t)e- '" ^ * " . 

For \x\ < \a~\t, we divide the analysis into cases, s £ [Q,t/2] and s £ [t/2,t]. For s E [Q,t/2], we observe the 
estimate 

(x-a-(t-»))2 

e s^T?^^) (1 + s)~^ 

(6.29) 



< C 



(x_a, (t-^.))^ 



i:^^ (l + s)"^+e {l + \x- ajtiy^ 



For the first estimate in H6.29|l . we proceed similarly as in H6.26(l . while for the second we have an estimate 

by 

Cit-\1 + \x- ajt])-^ [ e ds < C{1 + t)-^/^{l + \x - a-t\)-\ 

Jo 

which is sufficient for t > |a;|/|a^|. For s £ [t/2,t], we have 



(x-a {t-3)r (x-a, {t-3)r 

1/2^ M(t-s) <- r-i^i-i/z. 



(t_s)-V2g M(t-s) <C\x\~^'^e 2M(t-s) 
from which we obtain an estimate on H6.28|l by 



(x-a. (t-B)r 



C2|x|-i/2(l + t)-i / {t - s)-i/2e ^fc^ds < C|xri(l + 

Jt/2 



for which we have oberved the bound 



t (x-a-(t-3))2 

1/2, 



t/2 



(t - s)~i/^e 2M(t-=) < C. 
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The second estimate in (I6.27f) can be analyzed similarly. 

^8.1t^ i)). nonlinearity v(jj, s)ip{y, s). We next consider integration against the critical nonlinearity 
v{y,s)ip{y,s), which constituted the limiting estimate of the analysis in jl5|. Here, we refine the analy- 
sis of |15| both through refined estimates on v(y, s) and through the application of the approach of Liu 
|25[ described in Section 01 which takes advantage of improved decay for derivatives along the charac- 
teristic direction. It is precisely for this analysis that we must keep track of estimates on the characteristic 
derivatives [dt + ajdx)v. We proceed by dividing the integration over s as, 



[ G{x,t- s;y)(v{y,s)ip{y,s)) dyds 

J-oo ^ 'V 

/ G{x,t- s\y)[v{y,s)ip{y,s)\ dyds 
' ^ (6.30) 

/ G{x,t- s;y)[v{y,s)^{y,s)\ dyds 
■Jt J-oc ^ 'y 

G{x,t - s;y)(v{y,s)ip{y,s)j dyds. 



For the first integral on the right-hand side of Ht).30|) . we integrate by parts in y and use the supremum 
estimate ||t^(y, s)||loo < C(l -f s)~^/^ (valid for either of our estimates on v{y,s)) to obtain integrals of the 
form, 

{t-s)~^e STTt^T) (l + s)"^/'*e si^ dyds, 



10 J-oo 

for which we observe the equality 



M(t~s) g st: — — e Ml e M,(t-s)yy t > (6.31) 



(see Lemma 6 of |16)'). Integrating over y, we immediately obtain an estimate by 



(x-a, {t-s)- 



Ct-^l^ l^\t - S)-1/2(1 + g)-3/4e- "' '\n " ' ds. 







For \x — aj t\ < K^/t, any fixed K, we have kernel decay exp(— (.t — aj t)^/{Lt)) by boundedness, while for 



a; — a 



't\ > K^/t, with K sufficiently large, we have 



\x - a-{t - s) - a-s\ = \{x- ajt) + (aJ - a^)s\ > (1 - )\x - ajt], (6.32) 

and we again have kernel decay exp(— (a; — ajt)'^ / (Lt)). In either case, we obtain a final estimate by 

C(l + t)-3/4e- '^'ir' , 

which is sufficient by l|6.1|) . Similarly, for the third integral in l|6.3U|l . we obtain, upon integration in y 
precisely as above, an estimate by 

/* {t - s)-i/2(i + g)-^/v '^'°^ '^Mr^^'^ s. 

Proceeding similarly as above, we obtain an estimate in this case of the form 

(i + t)-3/V^^^^. 
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For the second integral in (|6.3U|) . we first consider the case j = k, for which, proceedingly similarly as above, 
we have integrals of the form 

Jvt 

which has been shown sufficient above. In the case j 7^ fc, the claimed estimate does not follow from such a 
direct method, and we employ a Liu-type cancellation estimate, as described in Sectional based on integration 
by parts in the characteristic direction. In order to clarify our analysis, we define the non-convecting variables 

9{x,t;y) = G^x,t,y- ajt) 
<f>{y,s) ^ (p'^.iy + a'^s,s) (6.33) 

where 

(ri-j,-o-(f-s))2 



G^(x,t;y)=rt-^/^e ' c ^ rj {1-^ / ^ AttPJ , (6.34) 

and (^^ is as in Ijl.f (We will consider corrections to at the end of the analysis.) In this notation, the 
second integral in H6.30|l becomes 

/ 9ix,t - s;y + aj{t - s))(v{y - a'^s,s)(f>{y - a^s,s)) dyds. (6.35) 
Setting ^ = y + a^(t — s), (|6.35|l becomes 



g(x, t - s; ^) {y{i - a - (t- s) - flj. s, s)(/)(^ - {t - s) - a^, s, s)^ ^d^ds. (6.36) 



Denoting by dr differentiation with respect to the second dependent variable of each function (effectively, a 
characteristic derivative dt + a^^^ on our original variables), we observe the differential relationship 

(^gix,t- s;^)V{S. ~ aj{t - s) - a^s,s)(/)(^ - aj{t~ s) - a~s,s)^^ 
= -9T{x.,t ~ s;^)V{£, - aj {t - s) - a]: s, s)(p{^ - aj {t ~ s) -a]:s,s) 

+ («7 ^ ^k)9ix^ t- s;^) (v{^ - aj{t- s)- a^s, s)<j){( - aj {t - s) - a^s, s)j^ ^'^^ 
+ 9{x, t- s;^) (v{£, -aj{t- s) ~ a~s, s)(j){S, - aj {t - s) - a~s, s)^ ^. 

Recalling that the case j = k has already been considered, we can rearrange ()6.37f) so that the integrand of 
(I6.36() can be written as 

g{x, i - s; (vi^ -aj{t- s)- a^s, s)(/)(^ -aj{t- s) - a^Ts, s)j ^ 

= («7 ^ (5(2^' * ~ 0^(C - aJ (t-s)- a^s, s)(/)(C - aj{t - s) - a^s, s)j ^ 

+ («7 - ak)~^9T{x,t- s;C)V{^ - aJ {t - s) - a~s,s)0(^ -aj{t- s)- a~s,s) 
- («7 ~ «fc)~^5(a;,i- s;C)(t^(C- a7(*~ ~ '^k^^^)^(^ - aJ {t - s) - a^s,s))^. 
For the integration over the first expression on the right-hand side of H6.38|l , we change the order of integra- 
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tion, and evaluate integration over s to obtain 

[g(x, t - s; ^)V{C -aj{t- s)- a^s, s)(j){^ ~ aj {t - s) - a~s, s) ) d^ds 
Vt J-oo ^ 

+00 

g{x, Vt; £,)y{^ -ajVi-a^{t-Vi),t- Vt)(t){£, - ajVt - a^{t ~ Vi),t - Vi)d£_ (6.39) 

-00 
+ 00 

g{x,t- Vt; - aj (t-V)- a~ V, Vt)<p{^ - aj [t-V]- a~ V, V)d£,. 

In each integral on the right hand side of 16.39fl . we employ the estimate 

\\V{-ML^<C{l + t)-^/\ (6.40) 
and proceed similarly as in |(n3U. For the first, we obtain an estimate by 

{VtT^''^e^^^ (1 + (i - Vt)y^''^e ml^) dS, 

-CO 

(x-a'7 Vt~a~ (f-\/t)2 

< ct-'/\i + ty'/'e '—^ , 

which gives an estimate by 

(see (|6.32() '). For the second integral on the right hand side of H6.39|l . we obtain an estimate by 

(t - Vr^^^e ^^<*-^) (1 + Vr^^'^e ^ d^ 

(a:-a-(t-Vt)-a-v/t)2 

< CT-i/2(i + vr'/'e — , 

which gives an estimate by 

(a;-a~t)^ 

C(l + t)-^/«e rf^. 



For integration over the second integrand on the right hand side of 16. 381 we have integrals 



Vt 



gr{x,t- s;^)V{C - a-j {t - s) - a^, s, s)(/>(^ - a - {t - s) - a^, s, s)d^ds 



<Ci / (i-s)-3/2e-MF^(H-s)-5/4e '—m d^ds (6.41) 

pt — y/t (x~-a~ {t — s) — a~s)-^ 

< C^t'^l^ \ {t - + s)-3/4e ' M. ' ds, 

where for the first inequality in (|6.41|) . we have employed the estimate (|6.40|l . while for the second we have 
used a calculation similar to (|6.31|) . In this last integral, we have three cases to consider, < < aJ , 
ttk < aJ < 0, and aJ < < 0, for which we focus on the second. (We recall that the case = aJ has 
already been considered above.) For |a;| > there is no cancellation between summands in (|6.6|l . and 

we obtain kernel decay, 

Ct-^/^{l + t)-^'^ ln(e + t)e-^^^^, 
while for |a;| < |a~|t, there is no cancellation between summands in (|6.8(l . and we obtain kernel decay 



Ci-i/2(l + t)-='/* ln(e + t) 



e ~* 
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In either case, the seeming blow-up as t — )■ can be ehminated by an alternative analysis in the case of t 
bounded. For \a~\t < \x\ < |a^|i, we divide the analysis into cases, s £ [^/t,t/2] and s £ [t/2,t— \/t]. For 
s £ [y/t,t/2], we observe through (|6.8|l the inequality 



(l + s)-'e-- 



< C 



- - 2 (6.42) 



(1 + sye — + (1 + la; - ajtir-'e 



where 7 = 3/4. For the first estimate in H6.42|l . we proceed as in the case < \aj \t, while for the second 
we obtain an estimate by 

<Ct-'/'-il + t)-'/'il + \x-aJt\)-^/\ 

which is sufficient. For s £ [t/2,t— ^/t], we observe through (|6.6(l the inequality with 7=1. For the 

second estimate in (|6.1(J|I . we proceed as in the case \x\ > \a'^\t, while for the first we have an estimate by 

C2t-'/^il+tr^/^\x-a-t\-^ e ' ds 



•'k 

< Cil + t)-^/^\x-a^t\-^ 



t/2 



which is sufficient for \x — t\ bounded away from 0. In the case of |a; — a^, t\ bounded, we proceed 
alternatively. 

The third term in can be regarded as the crucial piece, since it is here that we must keep track not 

only of estimates on v{y, s), but also on estimates of characteristic derivatives on v. We begin by expanding 
the characteristic derivative. 



g{x,t- s;^)(v{C-aJ{t - s) - s, s)(t){S, ~ aj {t ~ s) - a^s,s)j d^ds 

/ g{x,t-s;OV{^-aJ{t~s)-a^s,s)M^-aJ{t-s)-a^s,s)d^ds (6.43) 
ft—\[t /'+00 

/ g[x, t - s; ^)Vt{£, - aj{t- s) - a^s, s)(f>{(, - aj {t ~ s) - a^s, s)d^ds. 
'\/t J -00 

For the first integral on the right-hand side of H6.43|l . we employ the estimate H6.40II to obtain an estimate 

by 

/ {t - s)-i/2e-|T(A) (1 + + s)-3/2e ' M. d^ds 

V* J -CXI 

<Ct~'^' {l + s)-y'e ' " ds, 

which can be analyzed similarly as was Ht).41|l . For the second integral on the right-hand side of Ht).43fl . we 
note the following relation, also useful in calculations below, 

\Vr{^ - a^{t - s) - fljTs, s)| = K^f + flfc 9x)w|(C - aj [t - s), s) 

< Cs-\s'^\l + le - a-{t - s) - a^s\ + s'l^y^'^ 

+ ^ {\ + \i-a-{t~s)-a-s\ + s^'^)-^/^ 

af^O.l^k (6.44) 
+ Cs-^/\l + 1^ - a-{t s)\)-'/\l + le - a-[t ~ s)\ + s)-''^ 
X (1 + le - aT(t _ .)| + .i/^)-i/^/{,-,<(^_,- 

+ (1 + le - a7(t - s)|)-i(l + s)-i + (1 + le - a7(i - s)| -f s)-'/\ 
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For the first estimate in H(j.44|l . we use a suprcmum norm to obtain integrals 

t-Vt /•+00 (x-{)2 _ _ (J_a-(t-a)-a-a)^ 



{t - s)-^/^e~^ny^s^'^{l + s)-^e ■ — d^ds 

<Ct-'/' s-i/2(i + s)-ie ' ds, 

which can be analyzed similarly as was Ht).41|l . For the remaining estimates in (|6.44|) , the critical observation 
is that when integrated against the convecting diffusion kernel, 

e , 

they give increased decay in s, from which the claimed estimates can readily be observed. This completes 
the proof of jSUHKi)- 

Proof of (|3.18f ii)^. Integration against the nonlinearity s) has been considered in and we 
need only refine one estimate from that calculation in order to conclude our claim. The critical calculation 
regards integrals 

pt—'/t (^_„-(t_s)-a~s)2 

t-'/' it^sr\l + .s)-'/'e ' »^ " ds, 

Jt/2 

which was estimated in (15) by the expression, 

(l + t)-3/4(l + |:,_a^t|)-i/2. 

Here, we observe that alternatively, we may proceed by observing through (|6.6|l the inequality (|6.1Q() with 
7 = 1. For the first estimate in (|6.1()(l . we have an estimate by 

C2<-i/'(l+t)-i/2|a;-a-t|-i / e ' ' ds 

while for the second we have an estimate by 



t/2 



C2t-^^^il+t)-^e / {t^s)-^ds 



Jt/2 

< C(l + i)-3/2 in(e + ^)e- '^'^'" , 

either of which is sufficient. This completes the proof of H^.tHf ii)). 

Proof of (|3.18f iii^— (iv^'). x-derivatives. As will be clear from our analysis of characteristic derivatives 
just below, analysis of x-derivatives is almost precisely the same as the case of characteristic derivatives for 
which the direction of differentiation matches neither the convection rate of the kernel or the convection rate 
of the nonlinearity. Since we consider the case of characteristic derivatives in great detail, we will omit the 
case of x-derivatives. 

Proof of (IS^Stv)— (vi)), characteristic derivatives. We next develop estimates on characteristic 
derivatives {dt + a^dx), ^ 0, on the nonlinear interaction integrals. These estimates are the primary new 
contribution of the current analysis. 

Observing that G{x, 0; y) = dy{x)I by construction, we have, for af^ ^ 0, 

{dt+afdx)j J G{x,t- s;y)^<^{y,s) +T{ip,v,—S)y{y,s)^dyds 

= ^x,t)+T{^,v,^S)4x,t) (6.45) 



+ 00 ^ r 5y<5 

{dt + afdx)G{x, t- s;y) [^$(y, s) -f v, -^5)y{y-, s) 



dyds. 
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For C(i) > so that 

Hx,t)\ <cmM^:t) + M^,t)) 

\dMx, t)\ < Cit) [t-^^^M^, t) + ^2{x, t)) + ^^{x, t) + ^^{x, t) 
the first two terms on the right-hand side of (|6.45(l can be seen to be bounded by 
CE^{\ + t)-3/4 ^ (1 + 1^ _ a^t\ + ti/2)-3/2 



(6.46) 



t-\l + t)l/Vf * [X, t) + t-^''\4,{^ + XP2{X, t)) + ^3(2:, t) + X^i{x, t) 



Proof of (|3.18f vi^'). nonlinearity <!>(?;, s). The basic elements of our proof are more easily seen in 
the case of estimate (|3.18f vi')') — involving the diffusion wave nonlinearity s) — and so our approach will 
be to establish (|3.18r vi'l'l first and return to H^-lSf v')'). Proceeding from the right-hand side of H6.45|l . with 
nonlinearity s), we have integrals of the form 



t /'+00 



J~ao 



{dt + afdx)G{x, t- s; y)$(y, s)dyds. 



(6.47) 



In particuar, we focus on the nonlinearity {tp^ )j,, with iy9^ (j/, s) as defined in p.lO|l . and the leading order 
Green's kernel, 

&{x,t;y) = ct-^/'^e ^ , 



with c = {Ij Y'^ / yjAwPj (we will discuss corrections to G^{x,t;y) at the end of the analysis). We divide 
the analysis into three parts, 



t /'+00 



J-00 



{dt + fl; dx)G^ (x, t- s; v){iPk {y, s) )ydyds 



J -OQ 

t-Vt n+o 



(dt + ai d.j,)G^{x,t- s-y){ipi,{y,s) )ydyds 



(6.48) 



+ 



t 



{dt+ai dx)G\x,t- s; y) {<p^ (y, s) )ydyds 



t~\/i J- 



+00 



(9t-t-a, d.x)G^{x,t- s;y){ip^{y,s) )ydyds 



According to the estimates of Proposition [21 the case j = k does not occur (the diffusion waves have 
been chosen to eliminate precisely this case), leaving three cases to consider, l=jj^k, l~k^j, and 
I ¥^ 7^ k, j ^ k. We remark that it is in this third case that the characteristic derivative acts similarly as 
a derivative with respect to x only. 

Case 1: I — j ^ k. For the case I ^ j k, the characteristic derivative of G^ behaves like a time 
derivative of the heat kernel, and we have 

\{dt + a-dx)G\x, t~-s-y)\<C{t- s)-^/^e- °" . 

For the first estimate in H6.48|l . upon integration by parts in y, we have integrals of the form 



{t-s)-^e 



(j;-y-a (t-a)) 



' dyds, 



'Q J -00 

for which we observe the equality H6.31(l . Integration over y leads immediately to an estimate by 



' ds. 
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In the event that \x — Gj t\ < K\/i, we immediately have decay with scaUng exp{{x — Uj t)^/ (Lt)), while for 
|a; — ajt\ > K^/t we have 

\x - aj{t -s)- a^s\ = |(x - ajt) - (a^ - aj)s\ > (1 - °' )\x - ajt], (6.49) 

for which, with K taken sufficiently large, we again have decay with scaling exp((x — ajt)'^ / {Lt)). We have, 
then, an estimate on this term of the form 

Jo 

where in obtaining this last inequality we have reserved a small part of the kernel exp(— (x — aj (t — s) — 
s)"^ / {Mt)) for integration. Finally, the blow-up as t ^ Q can be reduced by putting derivatives on the 
nonlinearity for small time. For the third integral in (|6.48f) . we integrate the charactericstic derivative by 
parts to avoid blow-up near s — t. Observing the relation 



{dt + a - da:)G^{x,t- s;y) = -(9^ + a dy)G\x,t- s]y), 



we have 



t-Vt 



+°° . . 2 

{dt + aYdx)G^{x,t- s]y){ip'^{y,s) )ydyds 

-OO 

t /'+00 ^ _ ^ 

{-ds ~ a'^dy)G^{x,t~ s;y){(f^{y,s) )ydyds 



— OO 
+00 



CX2 



G^{x,Vi;y){^^{y,t~Vi)^)ydy (6.51) 
G^x,0;y){(pf, {y,t) )ydy 

+°° _ . 2 

G\x,t- s;y){ds + dy){ipl {y , s) )ydyds. 



For the first integral in H6.51|) . we have 

+ 00 (x-y-a^v/t) (y - a ~ (t - yT) ) ^ 

{Vt)-^/^e {l + {t- \/t))~3/2g M(t-vt) dy 

-OO 

< Gt-'^^{1 + {t- Vi))-'e , 

which gives an estimate by 

i-i/2(i + i)-ie-^^^, 

sufficient by an argument similar to H6.49|l . For the second integral in H6.51|l . &{x,0;y) is a delta function, 
over which integration yields an estimate by 

C(l + t)-3/2e" '°"A»"' . 

For the third integral in (|6.51() . we have, upon noting that the characteristic derivative is not along the 
direction of propagation of (/3j7 , integrals of the form 

/ {t~s)-^^^e W^^^ {l + s)-^/^e T^dyds 

t~\/t J — OO 

{x-a~ {t-s)-a'7 s)"^ 

<Ct-^'^ (l + s)-3/2e ' ds. 
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Proceeding similarly as in H().49|) and the surrounding estimates, we obtain an estimate by 

which has been shown sufficient in For the second integral in H6.48|l . we must proceed by taking 

advantage of increased decay for derivatives along characteristic directions. Recalling our definitions 
we can write this integral in the form 

I I grix^t- s\y + a~{t~ s)){(j}{y - a1s,sf)ydyds, 

which upon the substitution ^ ^ y + aj{t — s) becomes 

j*t—^/t /'+C30 

Proceeding similarly as in H().37|) . we write 

[gr[x, t - s; -aj{t- s) - a~s, sf^ ^ 
= -gTT{x,t~ s;^)0(^ -aj{t- s)- a^:s,s)^ 
+ (o-j - H )9T{x,t- s;O(0(C - aj{t-s) - 5,5)^)5 
+ grix, t - s; (^(^ -ajit- s)~ a^s, s)^)^, 
which can be rearranged as 

grix, t - s;^){(l){£, - aj {t ~ s) - a^s, sf)i_ 
= («7 - o-kT^ {9r{x,t ~ s]£,)4>{i - {t - s) - al s, sf^ 



(6.52) 



(6.53) 



where we recall again that the case j — k does not occur here. For integration over the first summand on 
the right-hand side of (|6.53|l , we exchange order of integration to obtain 



t—\ft /'+00 



+00 



00 
+ 00 



{gr{x, t - s; ^)(t){^ - a J {t - s) - a,, s, s)^^ ^d^ds 
grix, Vt- Ocj){£, -ajVt-a^it-Vt).t- Vifd^ (6.54) 
gr{x, t - Vt- C)(t){£. - a- {t^Vt)- a^Vt, Vtfd^. 



For the first integral in H6.54|l . proceeding similarly as in ()6.31|l . we estimate 

(Vt)"^/^e wr (1 + {t~ Vt)y^e M(t-v-t) 

-00 



< ct-i/2(Vt)-i(i + (i- Vi))-i 



/2. 



e 

for which we conclude in a manner similar to H6.32|) an estimate by 
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Proceeding similarly for the second integral in l(j.54|l . we estimate 



e ^ 



for which we conclude an estimate by 

For integration over the second summand on the right-hand side of (|6.53l) . we estimate 

t-Vt r+oo (x-e)^ (;^„-(t-.)-a^.)2 

/ (t-s) ^i^e (1 + s) dids 

\ft J —OQ 

< Ct-'^^ / it - + s)-i/2e ' ds. 

For this last integral we have three cases to consider, a^T < < aj, a^T < ^7 "-"^ 
which we focus on the second. In the case |a;| > |a^r|i, there is no cancellation between summands in (|6.6|l . 
and we have an estimate by 



t/2 



< Ct-^/^e 



while similarly in the case < \t, there is no cancellation between summands in H6.8(l . and we obtain 
an estimate by 

Ct-3/2e c^. 

For the case \cij\t < \x\ < \a^\t, we divide the analysis into cases s £ [\/t, t/2] and s e [t/2, i — Vt] (we take 
t large enough so that ^/t < i/2, proceeding alternatively for t bounded, in which we need not establish t 
decay). For s € [v^, t/2], we observe through H6.8(l the inequality H6.42|l with 7 = 1/2. For the first estimate 
in l|6.42|l . we proceed as in (|6.55|l . while for the second we have an estimate by 

Cii-5/2(i + 1^ _ aTt\)-'/^ f'^ 6- '°"°' < + i)-i/2(i + _ a7i|)-i/2. 

For s £ [t/2,t — \/t], we observe through (|6.6|l the inequality (|6.10|l with 7 = 3/2. For the second estimate 
in Hfi.l()|l . we proceed as in (|6.55|l . while for the first we obtain an estimate by 

C2t-'/^{l + tr^/^\x - a^t\-^/^ r^\t - g)-i/^e- '°"°' 'Mr'^^' rfg 



J t/2 

<Ct-^'\\^t)-^/^\x-alt\-^'\ 

which is better than the required estimate along the a^T directions (since j ^ k). For integration over the 
third summand on the right-hand side of (|6.53|l . we estimate 

t-Vt r+co (x-(f (;-a-(t-s)-a-a)2 

/ {t-sy^'^e (l-f s)"^e ~- d^ds 

Vt J -00 

< Ct-^'^ / {t - s)-\l + s)-3/2e ' " ds, 
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for which we can proceed almost exactly as in our analysis of the second summand on the right-hand side of 
(|6.53|) (the full analysis is omitted). 

Case 2: I = k ^ j . For the case I — k ^ j, the characteristic derivative of behaves like a time 
derivative of the heat kernel, and we have 

{dt + ald^)^^{y, s) = 0((1 + s)-'/')e-'-^^^. 

In general, our strategy for this case will be to integrate by parts, shifting the characteristic derivative onto 
ip. For the first integral in Ht).48|l . after integration by parts in y, we estimate 

{t - sy^/'^e mT^) (1 + sy^e dyds 



which is bounded by 

ct-^'\i + ty^/^e-'^^^^, 

sufficient for j ^ / by the argument of (|6.1|) . For the third integral in H6.48f) . we have 



,r"Vt (x-a~ (t-s)-ar s)^ 

<Ct-^'^ (t-s)-i(l + s)-i/2e ' ds, 



[dt + a^d^)G^ {x,t;y)dy{ip^ {y,sf)dyds 



{ds + aldy)G'\x,t;y)dy{Lpl{y,sf)dyds 

+00 

2 



GHx,0;y)dy{ip^iy,t)ydy (6.56) 



00 
+ 00 



+ / &ix,Vt;y)dy{cp^iy,t-Vt))^dy 

dy&{x,t- s;y){ds + a~dy){{ip'^f)dyds. 

t — \ft J —00 

For the first integral in H6.56|) . we have an estimate by 

[x, t)d,^- [x, t)\ < c{i + ty^'^-'-^^^, 

while for the second we estimate 

{\fty^''^e MVt (\j^{t~Vt)y^''^e M(t-viy' dy 

< cr^^^ii + it- Vt)y'^e , 

which gives an estimate by 

t-^/\i + ty^e — 

For the third integral on the right hand side of (|6.56() . we estimate 

(t - sy^e Mit-.) (1 _|_ s)-2g dyds 



t-Vt 



(x-a, (t-s)-a. s)^ 



< / {t - s)-i/2(I + 5)-3/2e ' ds, 

Jt-^/t 



(6.57) 



36 



Howard, Raoofi, and Zumbrun, Liu Estimates 



which gives an estimate by 



(x-o, ty 



For the second integral in (|6.48() . we first integrate by parts, moving the characteristic derivative onto (p^ , 
and then proceed by a Liu- type characteristic derivative estimate as described in Sectional We have 

t~\ft /'+C30 



t-Vt /.+■ 

k '-'V 



4-00 

&{x, Vt; y)dy{p\y. t ~ Vi)fdy (6.58) 

-OQ 

Gl{x,t-Vt-y)^\y,Vtfdy 

-OQ 

t-\/t !-+■ 
Vt 



G-'(x, t - s; y){ds + dy){ip''{y, s f)ydyds. 



For the first integral on the right hand side of (|6.58|) . we proceed exactly as in (|6.57|l to obtain an estimate 

by 

while for the second we estimate 

it - Viy^e Mi.t-'yt) (1 _|_ ^)-ie ^ dy 

-00 

which gives an estimate by 

An argument similar to that of H6.1|l shows that this last estimate is sufficient for j ^ I, which is the current 
setting. For the third integral on the right hand side of l|6.58ll . we proceed in terms of the non-convecting 
variables (|5.33f) . for which the integral can be re- written as 

t—^/t r+00 

/ 9ix,t- s]y + a~(t-- s)){dr{^iy - s, s)'^))ydyds, 
\ft J —00 

where we recall that dr represents differentiation with respect to the second argument of (j). Setting ^ = 
y + aj{t — s), this last integral becomes 

/ g{x,t- s]^)(dr {<!){£.- a~{t- s) - a~ 3,3)"^)) d^ds, 
Vt J -00 \ J 

for which we write 

(g{x, t - 3- Odrm - aj{t - s) - a^s, .s)^)) ^ 
= -9r{x,t- 3; ^)dr {(!){£. - aj{t - s) - a'^:s,sf) 
+ i^j - ak)g{x,t- s-£){dr{(j>{^ - aj{t - s) - a^s, sf))^ 
+ 9{x, t - s; i)drT{(l){^ ~aj{t~ 3) - a^3, s)^). 



(6.59) 
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For integration over the left hand side of (|6.59|) , we exchange the order of integration to obtain 

^+00 



5(2;, Vt; 09Ti(t>{C ~a-Vt-a-{t-Vt),t- Vtf)d^ 



+ 00 



g{x, t - Vi; OdrW^ -a-{t-Vi)- a^Vi, Vtf)d£, 



(6.60) 



For the first expression in (|6.6U|) . we estimate 



(\/i)"^/^e S5vr (1 + (< - \/i))"2g M(t-v^) 

-00 

< Ct-^'\\ + t)-3/2e ^-TiT , 



which gives an estimate by 

Ct-i/2(l+t)-3/2e 

For the second integral in H6.60|l . we estimate 

+ 00 ({-aT(t-Vt)-a-x/t)2 

(i-Vi)"^''^e «('-v^) (1 + Vi)"^e ^5vl 

-00 

< ct-^i\\ + ^ty^i^e ■' " — , 

which gives an estimate by 

sufficient for I ^ j- For integration over the first term on the right-hand side of H6.59|l . we estimate 



{t - s)"3/2g-|f(T^ ^ gy2^ -^ni — " rf^ds 

< Ct-^l^ \ (t ~ s)-\\ + s)-^l\ ' ds. 

We have three cases to consider here, < Q < aj , < aj < and aj < < 0, of which we focus on 
the second. In the case \x\ > |a^|t, we observe that there is no cancellation between summands in H6.6I) . and 
consequently we obtain an estimate by 

+ C,^-/2(l + 0-/^^-^ /"""(t - sr^ds (6.62) 

Jt/2 

< Ci-i/^(l + t)-ie 

Similarly, for |a::| < |aj|i, there is no cancellation between summands in (|6.8|l . and the same calculation gives 
an estimate by 

For \aj\t < \x\ < la^Tl^, we divide the analysis into cases s € [\/i,t/2] and s £ [t/2,t— \/t]. For s e [^/t,t/2], 
we observe through H6.8() the inequality H6.42|l with 7 = 3/2. For the first estimate in (|6.42() . we proceed as 
in 16.62|l . while for the second we obtain an estimate by 

Cit-3/2(i + 1^ _ a7i|)-3/2 e- """' "^»'°'° < Ct-'/^{l + t)-'/^il + \x- aTt\)-^/^, 
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which is sufficient. For s E [t/2,t — ^/i], we observe through the estimate (|6.1UI) with 7 = 1. For the 
second estimate in H6.1(J|I . we proceed as in (I6.62f) . while for the first we have an estimate by 

C2r^/^\x-a^t\-\l + t)-^^'' f °' ds < C{1 + ty'/^lx - a^t\~\ 

Jt/2 

which is sufficient since we are in the case t > \x\/\a^\. For integration over the third term on the right-hand 
side of Ht).59|l . we estimate 



Vt 



(6.63) 

rt-Vt (x-a7 (t-s)-or s)2 

< / (1 + s)-5/2e ' ds, 

which can be analyzed similarly as in the immediately preceeding case. 

Case 3: l^j,l^k, k^j. For the case I j, I k, k j, we proceed directly in the first and 
third integrals of (|6.48|) and through the method of Liu described in Section^ and elsewhere for the second 
integral in H6.48|l . In each case, we obtain precisely the reduced decay estimate for characteristic directions 
other than . This concludes the proof of Estimate ij^^Jvi)). 

Proof of (|3.18f v^^. nonlinearity T. The analysis of characteristic derivatives in the case of nonlin- 
earity J- constitutes the single most involved section of the paper. 

IS.l^ v)). Nonlinearity ip{y, s)v{y, s). The critical nonlinearity of is (/^(y, s)w(y, s), for which we consider 
integrals 

{dt + afd^)G^ {x,t - s;y)(Lp^{y,s)v{y,s))ydyds. (6.64) 

Case 1: I = j . For the case I — j, and for s G [0, t/2], we integrate by parts in y and employ a supremum 
norm on |?;| to arrive at integrals 

(t - s)~^e Mjt^) (1 + sY^I'^e st; dyds 
< / (t- s)-3/2(i + g)-3/v '°'"°^'Mr'''°°' ds. 



We have three cases to consider, a^T < < a^, a^T < < 0, and aj < a'^ < 0, of which we focus on 
the second. (We observe that for the nonlinearity (p{y, s)v{y, s), the case aj — arises.) For the case 
< a~ < 0, we first observe that in the case |a;| > |a^|i, there is no cancellation between summands in 
(|6.6|l . and proceeding similarly as in H6.62|l we determine an estimate by 

1/4 



Ci-"(l+t)^/*e 

In the case \x\ < \a~\t, there is no cancellation between summands in (|6.8|l and we obtain an estimate by 



C7r2(l+t)i/4e— i^. 

We note in particular that the case aJ — has been accomodated in this analysis. For the critical case 
lajli < \x\ < lofe 1^ (now j ^ k), we observe through l|6.8|l the inequality H6.42|l with 7 = 3/4. For the first 
estimate in H6.42II . we proceed similarly as in (|6.62|) . while for the second, we have an estimate by 

Cit-\l + \x- ajtl)-^/^ j ' e- °' °' ds < Ct-\l + t)-^/\l + \x - ajtiy''/^ 
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which is (precisely) sufficient for t > |a;|/|aj, |. For s e [t/2,t — ^/t], we do not integrate by parts in y, and 
consequently obtain integrals 

/ (i - s)-3/2e MTT^^^ s-i/2(^ ^ s)-5/4e ^f^dyds 

^ (6.66) 

Jtl2 



the last of which can be analyzed similarly as was (|6.65|) . We note that the form of the nonlinearity in 
(|6.66|) arises from a supremum norm on v (in the case that ip is differentiated) and from the observation 
that the estimates on Vy that do not decay at rate s~^/^ have spatial decay different from that of Lp so that 
when multiplied by (/s the combination decays at rate s"^/** or better. In the final case, s e — -v/i, t], the 
expression (t — s)^"^/^ is not integrable up to s = t, and we must proceed by integrating the characteristic 
derivative by parts. We have 

+00 

{dt +aYdj:)G-'{x,t~ s;y){ipl{y,s)v{y,s))ydyds 

t— \ft J —OQ 

/ {-ds - a^dy)G^ {x, t - s; y){ipl{y, s)v{y, s))ydyds 

t—y/t J —OQ 

&{x,0-y){ipl{y,t)v{y,t))ydy (6.67) 

-00 
+ 00 

G'{x, Vi;y){(p'k{y,t~ Vi)v{y,t- Vi))ydy 

-00 

t p-\-oo 

/ Gl{x,t- s;y){ds + dy){ipl{y, s)v{y, s))dyds. 

t — \/t J —OQ 

For the first integral on the right hand side of (I6.67|l . we have the immediate estimates 

\^i{x, t)v{x, t) I < cE^am + ty^'^e- ^^^^ 

\^'' (x, t)u, {x,t)\< CEomt-'^'il + t)-5/4e- ^^^^ , 

where the second of these follows from the observation above that the combination ip^{x,t)vx{x,t) decays 
at a rate faster than a product of the supremum norms. For the second integral on the right hand side of 
(|6.67|l . we estimate 

{Vty^^^e ^ (t - Vt)"^/^(1 + (i - %/t))~^/'*e "(t-v^) dy 

(x-a^Vt-aZ (t- v^))^ 

<Ct-i/2(i + i)-5/4e , 

which gives an estimate by 

Ct-i/2(l+i)-5/4e-^^^^if^. 

For the third integral on the right hand side of H6.67() . we have two nonlinearities to consider. We begin with 
integrals 

Glix, t - s- y)v(y, s)(ds + a~dy)if^{y, s)dyds, (6.68) 



for which we estimate 



/ {t-sy^e {l + sy'^'^e ~- dyds 

t— \Jt J ~OQ 

<Gt-'/^ f (t-s)-V^(l + g)-^/Si/^e- '^"°^'^"»"°'^°'' dg, 
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which gives an estimate by 
For the second nonhncarity, we have integrals 



3/2 (^-.*)^ 



Gl{x,t - s;y)(p''{y,s){ds + dy)v{y, s)dyds. (6.69) 



Observing again that the combination (f'^iy, s){ds + ajdy)v{y, s) decays at rate s^-'^(l + s)^'^/'*, we see that 
we can proceed as in the previous case. 

Case 2: I = k j. For the case I = k, we divide the analysis as in (|6.48|l into integrals 

t /' + 00 

/ {dt + a^dx)G^ {x, t - s; y){ipl {y, s)v{y, s))ydyds 

J-oo 

\/t /' + 00 

/ {dt + aJd^)G^{x,t-s;y){Lpl{y,t)v{y,s))ydyds 

t—yt r-\-oo 

/ {dt + aYd^)G^ (x, t - s; 7;)((^^ {y, t)v{y, s))ydyds 

t /'+C30 

/ {dt + a~d^)G^ {x, t - s; y){(p'^{y, s)v{y, s))ydyds. 
For the first integral in H6.7()() . upon integration by parts in y, we estimate 

/ (t-s)-3/2e mr^) {l + s)-^'^e ^ dyds 

J -oc 

JO 

which gives an estimate by 

sufficient for I ^ j. For the third integral in 16.70(1 . we do not integrate by parts in y, and consequently we 
can estimate 

{t-sy^e s~^''^{l + sy^'^e ~^ — dyds 

t— y/t J — OO 

< Gt-^l^ f {t - + sy'^^e- '"' ''mT^'^'^ s, 

which gives an estimate by 

ct-^^Hi + ty^e-'-^^^^, 

which is sufficient. For the second estimate in (|6.70|l . we integrate by parts as in (|6.67|) . moving the charac- 
teristic derivative onto the nonlinearity, and, when appropriate, moving the y derivative onto G^ . We have 
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{dt +aldx)&{x,t- s;y){ipl{y,s)v{y,s))ydyds 



t—\ft /'+00 

/ [-ds - a-^ dy)G^ {x,t - s;y)iip'^ {y, s)v{y, s))ydyds 

J-oo 
+00 

&{x,Vi;y){'Pkiy^t- ViHy,t- Vi))ydy (6.71) 



00 

-foo 



Gl{x,t-Vt; y)(p^ (y, Vt)v{y, \ft)dy 

t-Vt r+oo 



y 

—00 



Gl{x,t - s;y){ds + a,^dy){ip,^ (y, s)v{y, s))dyds 

For the first integral on the right hand side of (|6.71() , we proceed precisely as with the second term in (|6.67|l . 
For the second integral on the right hand side of (|6.71|) . we estimate 



< ct-\i + Vty^'/^e — — , 

which gives an estimate by 

suffificient for / ^ j. For the third integral on the right hand side of H6.71|) . we have two integrals to consider, 

+ / Gl{x,t- s;y)(p;:{y,s){ds + a^dy)v{y,s)dyds. 

J-Jt J -00 

For the first integral in H6.72(l . we employ the supremum norm on \v{y, s)\ and estimate 



/ Gl{x, t - s; y)v{y, s){ds + a^. dy)ip^ (y, s)dyds 

•J —00 



(6.72) 



Vt 



(t-sY^e [X + sY^l^e ~- — dyds 

< Ct-^l^ \ {t - s)-i/2(l + sY'l^e ' ds. 

JVt 



(6.73) 



We have three cases to consider for this last integral, a^, < < , a^, < < 0, and < < of which 
we focus on the second (we have already considered the case j = k). For \x\ > \a^\t, there is no cancellation 
between summands in H6.6|l . and we obtain an estimate by 

c{i + tY'^\-^^^^^, 

while for \x\ < \a~t\, there is no cancellation between summands in 1)6. 8|l and we similarly obtain an estimate 

by 

C(l + t)-"/8g-^^^^^i^^ 

sufficient for I ^ j. For the critical case \aj\t < \x\ < |a^|t, we divide the analysis into cases, s E [\/t,t/2] 
and s £ [t/2, t - y/i]. For s e [y/t^, t/2] we observe through the inequality with 7 = 3/2. For the 

first estimate in (|6.42|l . we proceed as above to obtain an estimate by 

11/8, 



C(l + t)-^V«e 
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while for the second we estimate 



Ct-\l + \x- a-t\)-^/^ / (1 + s)-i/4e ' Mt ds 

<C(l + i)-5/«(l + |.T-aT<|)-3/2^ 

sufficient for I ^ j. For s e [t/2,t — Vt], we observe through (|6.6|l the inequahty H6.1()|l with 7 = 1/2. For 
the second estimate in we immediately obtain an estimate by 

precisely as required, while for the first we estimate 



Ci-l/2(l + t)-7/4|^_„-i|-l/2 

<Cil + t)-y^\x-a^t\-'/^, 



e J"* ds 

t/2 



which is sufficient since t > |a;|/|aj, |. For the second integral in l|().72|l . we have five estimates on {{ds + aj^ dy)v\ 
to consider (see (j6.44ll V beginning with integrals 



Writing 



(x-a-a,- (t-s)) 



{t - s)-^e^^' 'Mi-s,^" + s)-l/4(l ^_ |y _ a-^l gl/2)-3/2g- ''''2°' 



a; 



- 2/ - a^- (t - s) = - {t-s)~ a^^ s) ~ {y ~ s), 



we observe the estimate 

(x-H-a7(t-s)) 



< c 



(x-H-a7(t-a))2 (x - g " ( t - a) - a " a )3 

(1 + I2; - a^sl + si/2)-3/2g- 



e Mit-s) g 

(x-a-a-ct-s))" 



(6.75) 



+ e (l + |x-a^- (t-s)-a^,s| + s^/^)" 



-3/2, 



(x-a, (t-s)-a. 



(x-a-a, (t-s))'' 



(x-a, (t-^.)-a. 



(t-s)"^e ' .s"^(l+s)"^/'*e J*^<*-=) (l + |y-afcS|+s^/2)"^/^e — dyds. 



For the first estimate in H6.75|l . we have integrals 

t — \[t /'+C30 

\ft J — 00 

(6.76) 

We have three cases to consider, < < aj , ^ ^^'^ '^J ^ "-k ^ 0, of which we focus 

on the second. For \x\ > \aj^\t, we observe that there is no cancellation between summands in (|6.6|l . and 
consequently that we obtain an estimate by 



(x-a. ty 



2 /.t/2 



s~^/^(l + s) " '^fc-^) ds 



(x-o-t)2 /-t-Vt 



(x-a, (t-.)-a.. 



t/2 



For |a;| < |aj|i:, we observe that there is no cancellation between summands in H6.6() . and consequently we 
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obtain an estimate by 

(x-a7t)2 l-t/2 (x-a-(t-s)-a-s)'' 

Cit-^e / e"' s-^l'^{l + sy^ds 

+ C2t-^il + t)-^e rf— / it-s)-^^^e~' '~^^T^> ds 

Jt/2 

which is sufficient for j ^ I. For \aj\t < \x\ < \a^\t, we observe the inequaHty 

g mF^I) 5-1(1 4. 5)-V2 



< c 



(x~a7t)2 (x-a (t-s)-aj^ 



g Z.t g M(t 

(x-a-(t-s)-a-.)2 



5-1(1 + (6.77) 



+ e "^^<'-=) |a;-aJ^|-i(l + |a;-a7t|)"i/2 
For integration over the first estimate in Ht).77|l . we have an estimate by 

(x-aTt)2 rt/2 (x-a7(t-a)-a-.)2 

Cit-^e / s-^^^{l + s)-^e^' ^^n^^ ds 

+ C2t-i(l + t)"ie ri— / (t-s)-i/^e"' ^^<^^^^ ds 

Jt/2 

< c{i + ty^^^e- '^'^'^ " , 

which is sufficient for I ^ j. For integration over the second estimate in (|6.77|) . we have an estimate by 
+ C2(l + t)-^'''\x - a-t\-\l + \x- a^t\)-^'^ f - s)-i/2e-''^ ds 



t/2 



< Ct-'/^\x - a-t\-\l + \x~ a-ir 



which, along with an alternative estimate in the case \x — ajt\ < is sufficient for I ^ j. We remark that 
the critical observation in this calculation was that since we require less decay along the j-characteristic, we 
use the same estimate (|6.77|l for both cases s € t/2] and s S [t/2,t — ^/i]. For the second estimate in 
H6.75fl . we have integrals 

t-Vi f+OO (x-»-a7 (t-s))2 

/ {t - s)-^e s~\l + s)-^/^{l + \x-a-{t-s)- a' s\ + si/2)-3/2 

Vi (6.78) 

X e Ms e Ms dyds. 

We have three cases to consider, < < aj , < ^ ^^^'^ '^J ^ ^ 0' of which we focus 
on the second. For \x\ > \a^\t, we observe that there is no cancellation between summands in H6.6|l . and 
consequently we obtain an estimate by 

l'i/2 (x — a~ (t — s) — a ~ s)'^ 

C,t-\1 + \x- a^t\)-'^' / s-i/2(l + s)-i/4e ' ms " ds 

+ C2t-\i + tr^'Hi + \^- «fc f~^\t - s)-'i^ds ^^-^^^ 



t/2 



<C7(l+t)-3/4(i + |a;_a-t|)-3/2^ 
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which is the required estimate since I = k. For the case \x\ < \t, we have no cancellation between 
summands in (|6.8|l and proceeding as in (|6.79|) we immediately obtain an estimate by 

C{l+t)-^^^{l + \x~aJt\)-^/\ 

For \a~\t < \x\ < |a^|t, we divide the analysis into cases s G [^/i,t/2] and s e [t/2,t— ^/i]. For s £ [^/i,t/2], 
we observe the estimate (|6.77|) . For integration over the first estimate in (|6.77|) . we have an estimate by 

Cit-'e / s-'/\l + s)-'e~' ' Ms ds 

< C{l+t)-^^^e 

while for integration over the second estimate in (|6.77|) . we have an estimate by 

Cit-^\x - ajt\-^{l + \x- ajt\)-^^^ / e ds 

Jt/2 

< C(l + t)-'/^\x - ajt\-\l + \x^ a-t\)~'/^ 

which, along with an alternative estimate in the case ja:— aji| < \/t, is sufficient for I ^ j. For s €E [t/2, t~\/i], 
we compute an estimate directly from H6.78|l 

C2t-\l + t)-^ {t-s)-^/^e ' Ms ds 

Jt/2 

<c{i+try\ 

which for t > \x\/\a'i^\ gives an estimate by 

C{l + \x\+t)-y\ 

For the remaining estimates in \ (ds + a^ dy)v\, we have decay with with a different scaling than in the diffusion 
wave. For example, we have terms of the form 

t-Vt I- + 00 ix-y-a-(t-s))^ {v-a7sf 

/ (t - s)-^e ^TTt^i 5-1/2(1 + + \y- a-s\ + s'l^Y^'^e ^dyds, 

where m ^ k. In such cases, we observe that for y near a" s we have exponential decay in s (which gives 
exponential decay in \/t), while for y away from fl-g, we have integrals of the form 



{t - s)-^e 5-1/2(1 + s)-i/2(l + s)-3/2e '-rf^dyds, 



which are better than previous cases (see (|6.73|) '). 

Case 3: I ^ k, I j . For the case I ^ k, I ^ j, we divide the analysis into precisely the same three terms 
as m if^Tnil . For the first integral in ^7U^, upon integration by parts in y, we estimate 

(i-s)-^/2g Mit-s) (l + s)-^/'*e ~- — dyds 







^0 

which is sufficient for I ^ j- For the third integral in (|6.7()|l . we estimate 

< (x-y-a-{t-s)f (B-a-s)2 

/ (t-s)-ie ^n^) 5-1/2(1 + s)-5/4g 

< C<-i/2 r (t - s)-i/2(l + g)-^/^e- '^"°^ "^^'''^'' rfa < C(l + i)-3/2e-^^^^, 
Jt-^/t 
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which is sufficient. For the second integral in H().7U|I . we begin with the case j — k, for which we write 

{dt +aYdx)G^ix,t- s;y)[ip''{y,s)v{y,s)) dyds 

Vt J-oo ^ 'V 

t/2 f+oo 

{dt+a'[d^)Gl{x,t- s;y)(p''{y,s)v{y,s)dyds (6.80) 



Vt 



{dt + aYdx)G\x,t- s;y)\^ip''{y, s)v{y, s)j ^dyds. 



It/2 

For the first integral on the right hand side of lt).80|l . we estimate 

(t - s)-^/^e S7(T^^^ (1 + 5)-5/4g if^ds 

^/t 

ft/2 (x-a-t)2 

< cr'/' / (< - s)-\i + s)-'^^s'/\ — ^ds 

Jvt 

< Ci(l + ^)~^/^e- '°' a°A" ds, 

which is sufficient for I ^ j- Similarly, for the second integral on the right hand side of (|6.8U|) . we estimate 

{t - s)-^e ^^^rr^ 5-1/2(1 + s)-^/^e ^^ds 

t/2 

J t/2 

< Ci(l + ty^^^e- ^'' ds, 

which again is sufficient for I ^ j. In the case j ^ k, we employ the non-convecting variables H6.33|l . along 
with 

gi{x,t;y + aj {t ~ s)) := {dt + aYdx)&{x,t- s;y), (6.81) 
for which we can write the second integral in (|6.7U|I as 

rt~\/t I'+oo 

/ / 9i{x,t- s;y + a~{t- s))((t){y ~ a';:s,s)V{y - a'^s,s)) dyds. 
Setting ^ = y + aJ {t — s), this becomes 

ft—\/t f+OO 

gi{x, t ~ s; O(0(^ -a~{t- s) - a^s, s)V{C -a~{t- s) - a^s, s)] dyds. (6.82) 

Vt J-oo \ J ' ' V 

Proceeding similarly as in 1)6.37(1 . we write 

[gi{x,i- s;0(t>{C - aj{t - s) - a',^s,s)V{^ - aj{t - s) - a',^s,s)j^ 

= -9ir{x,t- s;^)0(^ - aj{t - s) - a~s, s)V{C - aj{t - s) - a~s, s) 

+ («7 ~ a,'^)9i{x,t - s; ()(^(t){^ - aJ {t - s) - a'^: s , s)V {^ - aJ {t - s) - a^s,s)J ^ 

+ gi{x,t - s; ~ 0.J {t - s) - a'^: s , s)V {^ - aJ {t ~ s) - aj:s,s)^^. 

For integration over the left hand side of (|6.83|) , we exchange the order of integration to obtain 

gi{x, Vt; O'f'i^ -'a-Vi-a-{t-Vi),t- Vi)V{^ ~ aJ Vt - a' {t - Vi),t - V)d^ 



— oo 

+00 



gi{x, t - V; C)0(C -aJ{t-V)- a^V, V)V{^ aJ {t - V) - a^V, Vt)d£,. 



(6.84) 
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For the first integral in Ht).84l) . using the supremum norm of V, we estimate 

-00 

< ct-i/2(\/<)-i/2(i + (t _ vt))-^/Ht - Vty/^e — '—^ 

< C(f +t)~^/^e~ '^'^'" . 

For the second integral in Hti.84|l . again using the supremum norm of V, we estimate 

{t - V7)"^e"Af<*-^) (1 + Viy^^'^e M(t-v-t) 

-00 

(x-a~ (t--:/t)-a'r Vt) 

< Ct-^/\t ~ Vt)-i/'(l + Vt)-'>/\Vt)^'^e ' 

< C(l +t)-"/«e-^^T^, 

which is slightly better than required for I ^ j (we require t~^^^^). For integration over the first term on the 
right hand side of 16.83|l . we estimate 

{t - s)-2e- (1 + s)-5/4e M-s d^ds 

(6.85) 



00 

< r^\t - s)-3/2(i + g)-5/v/2e- '°'"''^"":;r°'°°'' dg. 

We have three cases to consider, a^T < < a^, a^T < *^ '^J ^ '^k which we focus on the 

second. For \x\ > la^Tl^, there is no cancellation between summands in H6.6|l . and we can estimate 

JVt Jt/2 (6.86) 

<C(l + t)-3/2e-^^^. 

In the case |a;| < |aj|<:, we have no cancellation between summands in l|6.8|l and similarly obtain an estimate 

by 

For the critical case \a~\t < \x\ < |a^|t, we divide the analysis into subcases s G [v^, t/2] and s £ [t/2, t—^/i]. 
For s S [\/t, t/2], we observe through (|6.8|) the inequality (|6.42|l with 7 — 3/4. For the first estimate in 16.42|l . 
we proceed similarly as in (|6.86|) . while for the second we estimate 

Ct-^l + \x- a-t|)-3/4 /* + ,s)-^/^.^/V '^'°'"^"'°'°°'' ds < Ci-3/2(l + \x~ a-t\)-^/\ 

which is sufficient for t > \x\/\a'^\. For s £ [t/2,t— \/t], we observe through l|6.6(l the inequality H6.I()|I with 
7 = 3/2. For the second estimate in (|6.10|l . we proceed similarly as in (|6.86|l . while for the first we estimate 



C2{l+t)-^'^\x~a-t\-^'^ f ' ^r"" ds < Cil + t)-'/^\x - a^t|-3/^ 



It/2 

For the third expression on the right hand side of (|6.83|) . we estimate 



(t-s)"^e Af(*-=)s"^(l + s)"^e ~= d^ds 

< / (t - s)-i/2s-i/2(i + g)-ie- '°'"°^'M:'"°''°' dg, 
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where we have once again observe the increased rate of time decay for the combination 

- aj(t - s) - fljTs, s)V{^ - aj{t- s) - a'^s, s)j ^. 

We have three cases to consider, < < aj , ^ ^ ^' ^^'^ '^J < '^k < 0, of which we focus on the 
second. For \x\ > \a'^\t, there is no cancellation between summands in Ht).6|l . and we can estimate 

Cit-^e ^ / s-^^^l + s)-^ds + C2t-\l+t)-^e ^ / {t~s)-^/^ds 

JVt Jt/2 (6.88) 

(x-a-t)2 

< C(l + i)-^/4e rf— , 

which is sufhcient for k ^ k. For \x\ < \aj\t, there is no cancellation between summands in l|6.8|l and we 
obtain an estimate by 

C(l + t)-^/4e 

For \aj\t < \x\ < \a'i^\t, we divide the analysis into cases, s £ [^/t,t/2] and s e [t/2,^/t]. For s e [\/t,t/2], 
we observe the inequality (|6.42|l with 7 = 3/2. For the first estimate in H6.42|l . we proceed as in Ht).88|) . while 
for the second we estimate 

J\/i 

which is precisely enough for I ^ j. For s e [t/2,t — \/t], we observe through (|6.6|l the inequality H6.10|l with 
7 — 1/2. For the second estimate in (|6.10|l . we proceed similarly as in (|6.88|l . while for the first we estimate 

C2t-\l + t)-^\x-a-t\-^^^ [ "^e" '^'"' ^'mT^''^ ds < C{l + t)-^/^\x~a^t\-^/^, 

Jt/2 

which is sufficient for I ^ k. This ends the proof of (|3.18r vi)) for the leading order convection kernels 
GHx,t;y). 

fS.lS^ v)). Nonlinearity {v{y^ sy)y. We next consider integrals 

/ {dt + afd^)& (x, t - s; y){v{y, sf)ydyds 
J -00 

/ {dt+a^d^)G^{x,t-s-y){v{y,s)%dyds (6.89) 

^-00 

t n-\-OQ 

/ {dt + afdx)G^{x,t- s;y){v{ij,sf )ydyds, 

t/2 J ~oo 

for which we have two cases to consider, I ^ j and I = j. 

Case 1: I ^ j. For the case I ^ j, we integrate the first integral on the right hand side of (|6.89|) by parts 
in y to obtain integrals 

rt/2 f + OD 

dy{dt + afdx)& {x,t - s\y)v{y,s)'^dyds. 

Taking supremum norm on one v{y, s), we have two cases two consider, one for each estimate on v. For the 
first, we have 

t/2 /.+OC (^_„_„--(t_s))2 

{t ~ s)-3/2g MTT^^ ^ s)-3/4(i + |y _ a-s\ + s^'^y^'^dyds, (6.90) 

where < 0. We observe through (|6.3|) the inequality 1)6. 4|l with 7 = 3/2. For the first estimate in 16. 4|) 
we have integrals 

/ / (i-s)-3/2e-' e il + s)-^/\l + \y-a^s\ + s^/^)-'^/^dyds. (6.91) 
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We have three cases to consider a^, < < s, a^, < < 0, and < a^, < 0, of which wc focus on the 



second. In the case s e [0, \/t], we have 

(x-a-(t-a)-a-a)2 



e M{t-s} < —t ^ 

from which we obtain an estimate by 

Jo 

sufficient for I ^ j. For s £ [^/t,t/2], we begin with the case |x| > \a^\t, for which there is no cancellation 
between summands in (|().6|l . and we obtain an estimate by 

Cit^^/^e — / {1 + sy^e ^^(^^^^0 ds. (6.92) 

JVi 

In the event that j = k, we obtain an estimate by 

C(l + t)-3/2 ln(l + t)e~ '°""" , 
which is sufficient for k = j ^ I, whereas in the event that j ^ k, we obtain an estimate by 

C(l + t)-3/2e 

Proceeding similarly in the case \x\ < \a~\t^ we obtain an estimate by 

C{1 + t)-3/2 ln(l + t)e- '^'^t" , 

which is sufficient for / ^ j. For \aj\t < \x\ < \a'i^\t, we observe through (|6.8|) the inequality Ijti.yf) with 
7 = 3/4. For the first estimate in (|6 9|) . we proceed as in Ht).92|l . while for the second we estimate 

rt/2 (x-a7(t-.)-,-a)2 

Cit-^/^{l + \x- ajt\)-^ / e ^^(^^^0 ds < Cr^l + \x ~ aJt\)-\ 
which is sufficient for t > \x\/\a'f^\ and I ^ j. For the second estimate in H6.4|l we have integrals 

t/2 n+OC (x-y-aT(t-s))2 

/ {t - s)-3/2e mi^) (1 + s)-3/4(i + |y _ a-s\ + \x- aj {t - s) - a- s\ + s^'^y'^'^dyds, 

for which we observe through H6.8|l the inequality H6.13|l with 7 = 3/4. Integrating the exponential kernel, 
for the first estimate in H6.13|l we obtain an estimate by 

Cit-\l + |x - a-t\)-^'^ f'\l + sr^/Us < Ct-^/\l + \x- a7^|)-=^/^ 

Jy/t 

while for the second estimate in (|6.13() . we obtain an estimate by 

/•t/2 

Cit-\l + \x- a-t|)~'/^ / {l + \x- a-t\ + s^/^)-^/^ds < Ct-\l + |x - aJtiyK 

J \/t 

For the second estimate on v{y, s), H6.9()(l is replaced by 

''^ f {t - s)-3/2e-*^""M(t'=7'" (1 + |y|)-i/2(i + |y| + s)-5/4(^ ^ |y| ^ s^/^yV^dyds, (6.93) 
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for which we observe through l(j.8ll the inequahty lt).17|l . For the first estimate in (|6.17|) . we have integrals 

{t ~ s)-3/2e-^ mr^) e (1 + |y|)"'/'(l + |y| + s)-^/*(l + \v\ + s^/^y^^^dyds, 

for which we have two cases to consider, aj < and aj > 0. For the case aj > 0, there is no cancehation 
between a; and o,J{t — s) and the claimed estimate can be deduced in straightforward fashion. For the case 
aJ < 0, we first consider the subcase |a;| > |aj|t, for which there is no cancellation between summands on 
the right hand side of 

X - aJ {t - s) = {x ~ ajt) + ajs, 
and we immediately obtain an estimate by 

Cit-^'^e — / {\ + s)-^ds <C{l+t)-^'^\n{e + t)e (6.94) 



which is sufficient for I ^ j. For < \t, we observe the inequality ()6.19f) . For the first estimate in H6.19|l . 
we proceed similarly as in (|6.94|) . while for the second we estimate 



Cit-3/2(i + _ a7i|)-3/2 + s)i/2e-'^ < C{1 + t)-^/^{l + \x - aJt\)-''^'\ 



which is sufficient for I ^ j. For the second estimate in (|6.17|l . we have integrals 

(t - s)-3/2e-'' V(i-V (1 + \y\ + 1^ _ ajit - + |y| + \x ~ aJ [t ~ s)\ + s)-5/4 

X (1 + |y| + \x- a-(t -s)\ + s^/^)-^/^dyds, 

for which we have two cases to consider, < and > 0, and as before we need focus only on the former. 
For a~ > and \x\ > \aj\t, there is no cancellation between x — aJt and a~s, and we immediately obtain 
an estimate by 

Cit-\l + \x - ajt\)-^/^ f^\l + \x-aj{t- s)|)-3/4rfs < ^7(1 + + \x ^ aJt\)-^/\ (6.95) 



For < |t, we observe the inequality 16.23|l . For the first estimate in (|6.23() . we proceed similarly as in 
()6.95fl . while for the second we obtain an estimate by 

Cii-i(l + \x- a-i|)'3/2 /"'^"(^ ^ 1^ „ _ _,)|)-i/2^g < + i)-i/2(i + _ a-i|)-3/2, 

Jo 

which is sufficient for I ^ j. For the case s e [t/2, t], we need to shift the characteristic derivative onto the 
nonlinearity. We accomplish this precisely as in H6.51|l . computing 

t /' + 00 

/ {dt + aYd^)&{x,t- s]y){y(y,sf)ydyds 

t /' + 00 

/ {-ds - a^dy)G\x,t- s;y){v{y,sf)ydyds 

t/2 J ~oc 

Gl{x,tl2-y)v{y,t/2fdy (6.96) 

-OO 

&{x,0;y){v{y,tf)ydy 

-OO 

t /'+00 

/ Gl{x, t - s; y){ds + aYdy)v{y, s^dyds, 

t/2 J~oc 
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where in the first and last integrals on the right hand side we have additionally integrated by parts in y. For 
the first integral on the right hand side of (|6.96|) . and for the first estimate on v{y, s), we have integrals 



(6.97) 



+ 00 (ai-H-o" (t/2))2 

t-^e 5TW2^ (1 + \y~ a-it/2)\ + {1/2^/^''^^! + {t/2))~^'^dy 

-00 

<Ct-l/2(l+t)-3/2, 

which is sufficient for t > c\x\, for any constant c > 0. In the case \x\ > t/c, we observe the decomposition 

X - y - ^a^ t = [x - -a^ t - -Ofe t)-{y- -a^ t), 
through which we observe the inequality 



(x-H-a (t/2)) 



< c 



^^^{l + \y-^a-t\ + it/2)'/'ry^ 

i^-V-aJ {t/nf (x-aJ(t/2)-a~(t/2)) ~ 



{l + \y-^a-t\ + {t/2r^^)-'/^ 



(6.98) 



+ e-^^^^^^A^(l + |y - la^tl + \x - aJ{t/2) - a^{t/2)\ + [t/2f/^)-^/\ 



In the event that |a::| > t/c, for c sufficiently small, we have exponential decay in both |a;| and t for the first 
estimate in H6.98|l . while for the second we have an estimate by 

c,t-^/^{\+tr^'\\ + \x\+t)-^/\ 

which again is sufficient. For the second estimate on v{y, s) we have integrals 

rie-'^^'-fV^r" (1 + \y\)-^/\l + \y\ + {t/2))-''\l + \y\ + (t/2)i/2)-i/2(i + {t/2))-^/^dy 

(6.99) 

which is sufficient for t > c\x\, for any constant c > 0. The case > t/c can be analyzed similarly as in 
the immediately preceding case. For the second integral on the right hand side of H6.96|l . we observe that 
(x, 0; y) is a delta function with mass at x = y, and consequently, we have an estimate by 

C\v{x, t)v^{x, t)\ < C(l + t)-3/4 [t-l/2(^^(2.^ + ^2(X, t)) + t) + t) 

which is sufficient. For the third integral on the right hand side of ()6.96fl . and for the first estimate on 
{ds + Oj/ dy)v{y, s), we have integrals 



+00 



(x-y-a. (t-a)) 



t/2 J -OQ 

<ct~^/^{i + t)-^/ 



[t - s)-^e W^^l s-i(l + s)-i/2(l + \y - ais)\ + s^^^y^/^dyds 

(t - sy^e-^'^ "mc^-V" s-i/2(i + s)-i/4(i + |y _ ^-5)1 + s'/^)-y^dyds. 



t/2 J -00 



(6.100) 



This last integral has already been considered in the analysis of (|6.2(l . and we obtain an estimate by 

Ct-^/^{l + t)-^'^^^{x,t). 

We proceed in precisely the same manner for the second and third estimates on [ds + aj dy)v{y , s), obtaining 
an estimate by 



C 



t-^/^i4>i (x, t) + ^2(x, t)) + ^/.3(x, t) + t) 
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For the fourth estimate on (dg + dy)v{y, s), we have integrals 

(^t /'+00 



(x-j,-a, (t-s)) 



{t - s)-^e (1 + sY'^l^il + \y\ + + \y\Y^dyda 



It/2 J -oo 

In this case, we observe the inequahty 
e (1 + 



< C 



ATF^T) e mF^T) (1 + + e S7(t^^) (1 _|_ |y| + |x - a7(t - s)|)"^ 



For the first estimate in H6.102() . we have integrals 

(x-y-a7 (t-s))2 (x-aTCt-s)) 



t/2 -'-laj \s 



M{t-s) g 



(6.101) 



(6.102) 



^^'<*-) (1 + s)-3/4(i + Ij^I + + \y\)-^dyds, (6.103) 



for which we have two cases to consider, Gj < and > 0, and we focus on the former. For < and 
additionally \x\ > \a~\t, we have no cancellation between x — ajt and ajs, and we obtain an estimate by 

C2{l + t)-^/Hn{e + t)e-^-^-^ f \t - s)-^ds + ^2(1 + t)-^^^ ln(e + t)e' " f {t-s)-^'^ds 

Jt/2 Jt-1 



< C{1 + t)-^/^[\iiie + t)]2e 



(6.104) 



which is sufficient. For |a;| < |a \t, we observe the inequality (|6.20|l with 7 = 1/2, for which we estimate 



(x-a . (t-s))^ 



C{1 + <)-^/* ln(e + t){l + \x\)-^/^ / {t - sy^/^e ds < C(l + t)"^/* ln(e + t)(l + \x\)-^/^, 

Jt/2 



which is sufficient for t > \x\/\aj \. For the second estimate in (I6.102|l . we have integrals 

(x-H-a-(t-s))2 



{t - sy^e ssTt^^o (1 + \y\ + \x- aj{t - s)|)"^(l + sy'^'^dyds 



t/2 J ~\a-^ \s 

< c{i + ty'/^ f [t - s)-i/2(i + \x-a^{t- s)\y^ds < c(i + tyy\ 

Jt/2 

which is sufficient for |a;| < Kt, some constant K. For \x\ > Kt, we estimate 

r-t rO (x-y-a" (t-s))2 



{t - s)-ie s^rfc^l (1 + \y\ + \x- aj{t - s)|)"^(l + sy'^^'^dyds 

t/2J-\a-\s 

< c{i + tyy\i + \x\ + ty^ [ {t- sy^'^ds < c{i + ty''/*{i + \x\ + ty\ 

Jt/2 

For the final estimate on {ds + afdyy{y, s), we have integrals 

y y (x-H-a7(t-s))2 

/ / {t- sy^e T^Tit^i (1 + sy^/'^y + \y\ + sy'"^dyds 

Jt/2 J-|a7|s 

< c(i + ty^'^ f {t - sy^'^ds < c{i + ty^, 

Jt/2 



(6.105) 



(6.106) 



which is sufficient for |a;| < Kt. Since \y\ is bounded by \a-^ \s, for |a;| > Kt and K sufficiently large, we have 
exponential decay in both |x| and t. 
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Case 2: I ~ j . For the case I — j, we have additional decay at rate (i — s)~^/^, from which we immediately 
recover the claimed estimates. 

yj.lt^ v)). Nonlinearity \^^6\'^. We next consider integrals 

t /' + 00 

{dt + afd,)Gl{x,t-s-y){\ + s)-^e-^\y\dyds, 

where in this case additional y-derivatives on the nonlinearity give no additional decay. 
Case 1: I ^ j . For the case I ^ j, we have integrals 

t-l /.+00 (x-H-aJ(t-s))2 

it - s)-3/2e (1 + sy^e-'^^y^dyds 

(6.107) 







t P + OO (x-y-a^- (t-a))2 

/ {t-s)-'^e {1 + sy^e-'^^y^dyds. 

t-l J ~oo 

In either case, we observe the inequality Ht).27|l . For the first estimate in (|6.27|l . we have, upon integration 
of e~''l^l , integrals 

ft-l (x-a7(t-s))3 

(6.108) 

ft (x-a^. (t-3))2 ^ ' 

+ / (t-s)~^/^e {1 + sy^ds. 

i~ < 0, we first observe th„„ _ 



it-i 

Focusing as in previous cases on the subcase aj < 0, we first observe that for > \a~\t, there is no 
cancellation between x — ajt and aJ s, and we consequently have an estimate by 



(x-a, ty 



Cr^/^e- ' " ' / {l + s)-^ds + C2{l + t)-^e — / {t - s)-^^^e^^^ ds 



*/2 (6.109) 



< Ct-^/^[ln{e + t)]e — , 

where in this last inequality we have observed that for s G [i/2,t], we have 

for rji > 0. For < \aj\t, we divide the analysis into cases, s G [0,t/2], s e [t/2,t— 1], and s E [t — l,t\. 
For s e [0, i/2], we observe the inequality (|6.29|) . For the second estimate in (|6.29|) . we obtain an estimate 

by 

rt/2 (x-aT(t-a))2 

Cii-3/2(l + \x- ajt\)-^ / e ds < Cr^l + \x ~ ajtl)-^, 

Jo 

which is sufficient for t > |a::|/|a~| and I ^ j- For the first estimate in H6.29|l . we proceed as in H6.109|l . For 
s e [i/2,t], we observe the inequality H6.20|l with 7 = 1, for which we have an estimate by 

C2{l + t)-^\x\-^ f (t-s)-i/2e"* 2M(!-a)" ds + C3{l + ty^\x\~^ f ds 

Jt/2 Jt-1 

<Cil+t)-^\x\-\ 

Case 2: I = j . For the case I — j, we have additional decay at rate (i — s)^^/^, from which we immediately 
recover the claimed estimates. 

This ends the proof of estimate p.lSf v)) for the leading order convection kernel G^{x,t;y). 

l^'S.lt^ v)-(vi)). remainder estimates. In our proofs of H3.18r v)-fvi'l'l. we have considered only the leading 
order convection kernel 



&{x,t-y) = ct-^'^- ^ 



(x-y-a-t)^ 

e 
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We must also consider the remaining three scattering estimates (terms in S{x,t;y)) and additionaUy the 
remainder estimates R(x,t]y). Beginning with the remainder estimate, we focus our attention on the non- 
Unearity {(p''{y, s)'^)y (analysis of the remaining nonlinearities is similar). We have the decomposition 

t /'+00 

/ {dt + a-d,,)R{x, t - s; y){ip''{y, s)\dyds 

J-oo 

{dt + aYdx)R{x,t-' s;y){ip''{y,s)'^)ydyds 



J -oo 

t—^/t /' + 00 



(6.110) 



{dt + a; dx)R{x,t- s\y){(p^{y,sY)ydyds 



OO 



t p-\-OQ 

[dt + aYdx)R{x,t - s;y){ip''{y,sf)ydyds, 



-ty/Mt 



(6.111) 



r 

where for y < and aj < 0, 

,/ 

{dt + a-d.,)R{x, t;y)^Yl 0(e-''*)5,_a.t(-j/) + 0(e-''(l"-^l+*)) 

+ O ((t + l)-3/2e-^-+ + e-"!"!) r^t + Xfl^e-^--y—^tflMt 
+ ^ O ((i + l)-ie-"^^ + e-"l^l) r^{t + l)i/2e-(:'=-y-a, 

+ X{|,-t|>|,|}0((t + l)-l/2i-l)e-(.-ar(t^|yKT|))VMtg-r,.+ 

+ E X{|.-,>M}0((t + l)-V2i-i)e-(-4(*-I.K-|))VM*,-..-_ 

a^>0, a+>0 

In each estimate with decay t^'^l'^ or we have time decay better than that of characteristic derivatives of 
&{x, t] y), and we can proceed as in the above analyses. For terms J2j=i 0('3~''*)^a;-a*t(— y), we proceed as 
in the proof of l@J, in which the interactions arising from our relaxation of strict parabolicity are considered. 
The only genuinely new term is the exponentially decaying contribution, which has reduced decay in t. 
Focusing on this term, we integrate by parts in y, observing that in the Lax and overcompressive cases 
differentiation with respect to y improves t decay by a factor t"-'^/^ (this is a fundamental point of difference 
between the Lax and overcompressive cases considered here, and the undercompressive case). We have 
integrals 

e""!^! / / {t-s)-^e W^^^ {l + s)-^e ^r^dyds 



< Ce-^l-lt-i / (1 + s)-h'^'e ' ds 

Jo (6.112) 

Jo 



(x-a, 



< C(l +t)-3/''e-''l^le 



for which we observe that 



e 



(x-a . t)^ 
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For the third integral in H().110|) . we estimate 



e-"l^l / / (t-s)-i/2e (l + s)-3/2e dyds 



t-^ J- 



(x-a-(t-s)-a. sy 



< Ce-'n\-k-'/^ f (1 + g)-3/2gi/2g- '-- ds 

< Ce-''l^li-i/2(l + i)-ie / e""^ ^t?^ ds 



t-Vt 



< C(l + t)~^e~^'l^le~ '" , 



which is sufficient, precisely as above. For the second integral in H6.110|) . and for t large enough so that 
\/t < t/2, we estimate (integrating by parts in y ior s G [^/t,t/2]) 



e 



-"1^1 / / {t~s)-^e s?(^^^o {l + s)-^e ^Ti—dyds 

J\/i J -oo 

+ e-''l^l / / (<-s)-l/2g MTI^i^ (H_s)-3/2g jil—dyds 

Jt/2 J-oc 

For s e [Vi,t/2], we first observe that for the case j = fc, we have an estimate by 

which is sufficient, as observed above. For the case j ^ k, we recall the inequality (|6.42ll with 7 = 1/2. For 
the first estimate in (|6.42|l . we proceed as in l|6.114|l . while for the second we estimate 



Observing the estimate 



Cie-''l"lt-i(l + \x- aTt\)-^/^ f'^ e- '^'°' ' Mr'" °' ds 

< C{1 + i)-l/2g-r,|x|(^ ^ 1^ _ a-t|)-l/2. 



we observe that this estimate is sufficient. For s e [t/2, t — ^/t], we first observe that in the event that j = k, 
this last integral in H6.113|l provides an estimate by 

while in the event that j =/= k, we have an estimate by 

C(l + i)^^e-''l^l, 

either of which is sufficient. 

I^S.lS\f v}~fvi)). full scattering estimates. In the scattering estimates S{x,t;y) of Proposition |21 we have 
three corrections to G^{x,t;y), respectively from convection, reffection, and transmission. For convection, 
the term is 
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aj > 0, which satisfies 



\{dt + a-d,)Giix,t;y)\<C 



The estimate with t^"^/^ decay is precisely as in the case of , and can be analyzed similarly. For the 
estimate with exponential decay in |a;|, we can proceed precisely as with the exponentially decaying term 
arising in the analysis of R{x,t;y). The remaining corrections, from the reflection and transmission terms 
in S{x, t; y) can be analyzed similarly. 

In our analysis of G^{x,t]y), we employed the relation 

dxG^ {x, t; y) = -dyG^ {x, t; y) 

(see the argument of (|6.51|) . in which the characteristic derivative {dt + aj dx)& {x,t — s;y) is converted 
into a characteristic derivative in the variables of integration —{ds + aJ dy)G^ {x,t — s]y)). Designating the 
scattering term arising from reflection as 

Gfix,t;y)^cnt-'^^e {^^)' 
we observe that the analogous estimate is 



dxG%'{x,t;y) ^ -^dy&^\x,t-y) + 0{t-'/'e-'^^^^)e 



from which we see that the conversion from a characteristic derivative in variables x and i to a characteristic 
derivative in variables y and s takes the form 

{dt + ajdx)&^{x,t- s]y) = -{ds + a^dy)&^'' {x,t - s]y) 

_ (x-(a-/a-)y-a- (t-s))^ 

In this way, we again have precisely the improved decay required by the analysis, namely 

{ds + aldy)&i{x, t-s-y) = 0{{t - s)-3/2)e ^^J^^ . 

The argument involving transmission terms is entirely similar. 

Finally, in each case in which a characteristic derivative in x and t is shifted to one in y and s an 
exponentially decaying error term arises, such as (from reflection) 



0(t-i/2e-''l^l)e' 



In all cases, these terms can be analyzed as were the the exponentially decaying corrections in {dt + 
a-dx)R{x,t;y) (see (lOT^V 

This completes the proof of Lemma 13 □ 

7 Nonlinear integral Estimates II 

Finally, we complete the paper by establishing the nonlinear integral estimates of Proposition 01 
Proof of Lemma |4l To show (I 19(1 (i) we need to estimate 



7^*(a;)0(e-''«(*-^)),5,_a.(^„,)(-2/)/:f(y)T(2/,s)dyd, 



< 

(7.1) 

G I {e-''»<-'-''^)\T{-x + a*{t- s),s)\ds 







56 



Howard, Raoofi, and Zumbrun, Liu Estimates 



for the various sources T{y, s) arising in the bounds for J^, <I>. 
For example, the typical term 

\v\\vyy\{y,s) < + +7^2)(y,s) 

arising in the bounds for J- leads to sources 

Ti(y,s) = (l + s)-i/2V;i(y,s) 

^(l + ,)-l/2(l + |j,_„±i|)-3/2 

and 

< (1 + .)-V2(l + |y|)-l/2(i + |y| + ,)-l/2(i + ^ ,l/2)-l/2^ 

Substituting Ti in (|7.1|l . we obtain 











g-^„(*-s)(^ + 5)-^ (1 + - a,-i - (a* - 4){t - s)\)-^ds, 
which , by (|5.2() . is smaller than 

-voit^^)^l + + 1^ _ ari|)-i(|l + (a* - ar)(t - s)\f^ds, 



t 

e 



(7.2) 



(7.3) 



(7.4) 



(7.5) 



which in turn (absorbing t ~ s powers in e "2 ' ) is smaller than 

(1 + |x - art\)-i f e-^(*-^)(l + ^Y^ds < C{1 + \x - art\)-i{l + t)-^'^ 
Jo 

Similarly, substituting T2 in (|7.1|l . observing that 

(1 + \y\ + <V2)-l/2 _ + |y|)-l/2^ (1 ^ ^)-l/4}^ 

and following the same procedure, we obtain an estimate of 
C(l + t)-i/2(l + \x\)-^/\l + \x\ + t)-i/2 min{(l + \x\)-^/^, (1 + t)'^/^} ^ 

C(l + t)-'/\l + \x\r'/'{l + \X\ + t)-'/\l + \X\ + il/2)-l/2 

(7.6) 

where, in the final step, we have used the fact that, for x{x,t) — 0, 

(1 + |x|)-l/2(l + \x\ + <)-l/2(l + \x\ + <l/2)-l/2 _ (1 + ^ ^)-3/2 < C'V'l(x,t). 

Bounds for other cases follow similarly. □ 
Acknowledgements. The authors were partially supported by the National Science Foundation, under 
grants DMS-0500988 (Howard), and grants DMS-0070765 and DMS-0300487 (Raoofi and Zumbrun). 

References 

[1] J. Alexander, R. Gardner, and C. K. R. T. Jones, A topological invariant arising in the analysis of 
travehng waves, J. Reine Angew Math. 410 (1990) 167-212. 



57 



Howard, Raoofi, and Zumbrun, Liu Estimates 



L. Q. Brin, Numerical testing of the stability of viscous shock waves, Ph.D. dissertation, Indiana Uni- 
versity, May 1998. 

L. Q. Brin, Numerical testing of the stability of viscous shock waves, Math. Comp. 70 (2001) 235, 
1071-1088. 

L. Brin and K. Zumbrun, Analytically varying eigenvectors and the stability of viscous shock waves, to 
appear. Mat. Contemp. (2003). 

T. Bridges, G. Derks, and G. Gottwald, Stability and instability of solitary waves of the fifth-order KdV 
equation: a numerical framework, Phys. D 172 (2002), no. 1-4, 190-216. 

J. W. Evans, Nerve Axon Equations I IV, Indiana U. Math. J. 21 (1972) 877-885; 22 (1972) 75-90; 22 
(1972) 577-594; 24 (1975) 1169-1190. 

H. Frcistiihler, Some results on the stability of non-classical shock waves, J. Partial Diff. Eqs. 11 (1998) 

23-38. 

H. Freistiihler and P. Szmolyan, Spectral stability of small shock waves. Arch. Ration. Mech. Anal. 164 
(2002) 287-309. 

H. Freistiihler and K. Zumbrun, Examples of unstable viscous shock waves, unpublished note, Institut 
fiir Mathematik, RWTH Aachen, February 1998. 

R. Gardner and K. Zumbrun, The Gap Lemma and geometric criteria for instability of viscous shock 
profiles. Comm. Pure Appl. Math. 51 (1998), no. 7, 797-855. 

J. Goodman, Nonlinear asymptotic stability of viscous shock profiles for conservation laws. Arch. Ra- 
tional Mcch. Anal. 95 (1986), no. 4, 325-344. 

P. Howard, Pointwise methods for stability of a scalar conservation law. Doctoral thesis (1998). 

P. Howard, Pointwise estimates on the Green's function for a scalar linear convection- diffusion equation, 
J. DifTcrcntial Equations 155 (1999) 327 36. 

P. Howard, Pointwise Green's function approach to stability for scalar conservation laws. Comm. Pure 
Appl. Math. 52 (1999) 1295-1313. 

P. Howard and M. Raoofi, Pointwise asymptotic behavior of perturbed viscous shock profiles. Preprint 
2005. Available: www.math.tamu.edu/^phoward/mathpubs.html. 

P. Howard and K. Zumbrun, Stability of undercompressive shock profiles, to appear, J. Differential 
Equations. 

J. Humphcrys and K. Zumbrun, An efficient shooting algorithm for Evans function calculations in large 
systems, Preprint 2005. 

J. Humpherys and K. Zumbrun, Spectral stability of small amplitude shock profiles for dissipative sym- 
metric hyperbolic-parabolic systems, Z. Angew. Math. Phys. 53 (2002) 20-34. 

J. Humpherys, B. Sandstede, and K. Zumbrun, Efficient computation of analytic bases in Evans function 
analysis of large systems, Preprint 2005. 

C. K. R. T. Jones, Stability of the traveling wave solution of the FitzHugh-Nagumo system. Trans. 
Amer. Math. Soc. 286 (1984), no. 2, 431-469. 

T. Kapitula and B. Sandstede, Stability of bright solitary-wave solutions to perturbed nonlinear 
Schrodinger equations, Physica D 124 (1998) 58-103. 

S. Kawashima, Systems of a hyperbolic-parabolic com,posite type, with applications to the equations of 
magnetohydrodynamics, thesis, Kyoto University (1983). 



58 



Howard, Raoofi, and Zumbrun, Liu Estimates 



[23] S. Kawashima and A. Matsumura, Asymptotic stability of traveling wave solutions of systems for one- 
dimensional gas motion, Comm. Math. Phys. 101 (1985), no. 1, 97-127. 

[24] S. Kawashima, A. Matsumura, and K. Nishihara, Asymptotic behavior of solutions for the equations of 
a viscous heat- conductive gas, Proc. Japan Acad. Scr. A Math. Sci. 62 (1986), no. 7, 249-252. 

[25] T.-P. Liu, NonUnear stabiUty of shock waves for viscous conservation laws. Memoirs AMS 56 (1985), 
no. 328. 

[26] T.-P. Liu, Interaction of nonlinear hyperbolic waves, in: Nonlinear Analysis, Eds. F.-C. Liu and T.-P. 
Liu, World Scientific, 1991 171-184. 

[27] T.-P. Liu, Pointwisc convergence to shock waves for viscous conservation laws, Comm. Pure Appl. Math. 
50 (1997), no. 11, 1113-1182. 

[28] A. Majda, The stability of multi-dimensional shock fronts: a new problem for linear hyperbolic equa- 
tions, Mem. AMS 275 (1983). 

[29] A. Majda, The existence of multi-dimensional shock fronts, Mem. AMS 281 (1983). 

[30] A. Majda, Compressible fluid How and systems of conservation laws in several space dimensions. 
Springer- Verlag, New York (1984) viii-|- 159 pp. 

[31] A. Majda and R. Pego, Stable viscosity matrices for systems of conservation laws, J. Diff. Eqs. 56 (1985) 
229-262. 

[32] C. Mascia and K. Zumbrun, Pointwise Green's function bounds and stability of relaxation shocks, 
Indiana Univ. Math. J. 51 (2002), no. 4, 773-904. 

[33] C. Mascia and K. Zumbrun, Pointwise Green function bounds for shock profiles with degenerate viscosity. 
Arch. Rat. Mech. Anal. 169 (2003) 177-263. 

[34] C. Mascia and K. Zumbrun, Stability of large-amplitude 

[35] A. Matsumura and K. Nishihara, On the stability of travelling wave solutions of a one- dimensional 
model system for compressible viscous gas, Japan J. Appl. Math. 2 (1985), no. 1, 17-25. 

[36] R. Plaza and K. Zumbrun, An Evans function approach to spectral stability of small- amplitude viscous 
shock profiles, preprint (2002). 

[37] M. Raoofi, LP asymptotic behavior of perturbed viscous shock profiles, to appear, J. Hyperbolic Differ- 
ential Equations. 

[38] M. Slemrod, Dynamic phase transitions in a van der Waals fluid, J. Differential Equations 52 (1984), 
no. 1, 1-23. 

[39] K. Zumbrun, Dynamic stability of phase transitions in the p-system with viscosity-capillarity, SIAM J. 
Appl. Math. 60 (2000) 1913-1924. 

[40] K. Zumbrun, Stability of large amplitude shock waves of the compressible Navier-Stokes equation. 
Handbook of Fluid Dynamics, Volume IV (2004). 

[41] K. Zumbrun, Refined Wave-tracking and Nonlinear Stability of Viscous Lax Shocks, Methods Appl. 
Anal. 7 (2000) 747-768. 

[42] K. Zumbrun and P. Howard, Pointwise semigroup methods and stability of viscous shock waves, Indiana 
Univ. Math. J. 47 (1998) 741-871. See also Errata, Indiana Univ. Math. J. 51 (2002) 1017-1021. 

[43] D. Serre and K. Zumbrun, Viscous and inviscid stability of multidimensional planar shock fronts, Indiana 
Univ. Math. J. 48 (1999), no. 3, 937-992. 



59 



Howard, Raoofi, and Zumbrun, Liu Estimates 



Peter HOWARD 
Department of Mathematics 
Texas A&M University 
College Station, TX 77843 
phoward@math.tamu.edu 
and 

Mohammadrcza RAOOFI 

Max Planck Institute for Mathematics in the Sciences 

Inselstrae 22-26 

D-04103 Leipzig, Germany 

mraoofi@indiana.edu 

and 

Kevin ZUMBRUN 
Department of Mathematics 
Indiana University 
Bloomington, IN 47405-4301 
kzumbrun@indiana.edu 



60 



